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Abstract 

To adopt a practical method to calculate the action of geometrical operators on quantum states is a crucial task in loop 
quantum gravity. In the series of papers, we will introduce a graphical method, developed by Yutsis and Brink, to loop quantum 
gravity along the line of previous works. The graphical method provides a very powerful technique for simplifying complicated 
calculations. In this first paper, the closed formula of volume operator is derived via the graphical method. By employing suitable 
and non-ambiguous graphs to represent the acting of operators as well as the spin network states, we use the simple rules of 
transforming graphs to yield the resulting formula. Comparing with the complicated algebraic derivation in some literatures, our 
procedure is more concise, intuitive and visual. The resulting matrix elements of volume operator is compact and uniform, fitting 
for both gauge-invariant and gauge-variant spin network states. 

PACS numbers: 04.60.Pp, 04.60.Ds 


1 Introduction 

As a non-perturbative approach to quantum gravity, loop quantum gravity (LQG) has made considerable achievements (see 
[1,2] for review arcticles, and [3, 4] for books). This theory rigorously enforces the lesson of general relativity and is built 
on a strict mathematical fundament. In LQG, the quantum kinematical Hilbert space was successfully constructed with the spin 
network states as its orthonormal basis. The elementary operators are the holonomy and flux operators. By suitable regularization 
schemes, quantum geometric operators, such as the length, area, and volume operators corresponding to their classical quantities, 
were well defined on the kinematical Hilbert space 'LQm [5, 6, 7, 8, 9, 10, 11]. The area operator plays an important role in 
computing the entropy of black hole. The volume operator is a cornerstone on which some physical interesting operators, for 
instance, the Hamiltonian constraint operator determining the quantum dynamics of LQG can be constructed. 

There are two versions of volume operator in the literature. The first one, based on ‘external’ regularization scheme, was 
introduced by Rovelli and Smollin in loop representation [8], and re-obtained in the connection representation [10]. The second 
one, based on the ‘internal’ regularization, was firstly defined by Ashtekar and Lewandowski [10]. In [11], Thiemann presented 
a rather short and compact regularization procedure to re-derive the second version of volume operator. Playing a crucial role 
in LQG, the spectra of the volume operator is pursued. Certain matrix elements of volume operator were calculated in the 
framework of loop representation by using graphical tangle-theoretic Temperley-Lieb formulation in [12], Then they were also 
derived in connection representation by a rigorous but tedious algebraic method in [1 1, 13], and their special case was re-derived 
using generalized Wigner-Echart theory [14]. Although those components of the volume operator are rigorously defined, the 
computation of their actions on the spin network states are difficult. The main reason is the following. The volume element 
operator at a vertex v of a graph y reads V v = ^/\q v \. Although the matrix elements of q v can be calculated using recoupling 
theory, the matrix has no obvious symmetries and hence is difficult to be diagonalized analytically for the case that the dimension 
of the matrix is bigger than nine. On the one hand, the derivation of closed formula in [13] is rigorous. But there is no universal 
formula with so tedious and abstract method. Hence, it is desirable to develop certain rigorous, but concise and intuitive method 
to calculate the matrix elements of volume operator, as well as other geometric operators. 

Graphic calculus has been introduced in LQG in a few papers (see e.g., [12, 15,16,17,18, 19]). But a concise, non-ambiguous 
and complete scheme is still desirable. In this first paper of the series, we will introduce a graphical method to LQG, which is 
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based on the Brink’s graphical method [20] and its suitable extension 1 . This method consists of two ingredients, graphical rep¬ 
resentation and graphical calculation. We will represent the algebraic formula by corresponding graphical formula in an unique 
and unambiguous way. Then the graphical calculation will be performed following the simple rules of transforming graphs, 
corresponding uniquely to the algebraic manipulation of the formula. The basic components of the graphical representation and 
the simple rules of transforming graphs are presented in Appendix A. A central goal of this paper is to derivate the closed formula 
for the matrix element of the volume operator, which involves only the flux operator, based on our rigorous graphical method. 
Comparing to the algebraic method, our derivation is obviously more compact and simple. Our analysis shows that the formula 
of the matrix elements for certain case in [13] is also valid for other cases and hence can be regarded as a general expression. 
In the second paper of the series [22] , we will consider the actions of the gravitational Hamiltonian constraint operator and the 
inverse volume operator, depending also on the honolomies, on spin network states in the graphical method. In order to obtain 
the on-shell anomaly-free quantum constraint algebra, one has to employ degenerate triangulation at the co-planar vertices of 
spin networks in the regularization procedure of Thiemann’s Hamiltonian. This problem has been overcome by a new proposed 
Hamiltonian constrain operator in [23]. 

This paper organized as follows. Section 2 is devoted to review briefly the elements of LQG. In section 3, we will introduce 
the graphical method to LQG. In section 4, we will derive the closed formula for the matrix element of the volume operator by 
the graphical method. It is shown how the simple rules of transforming graphs tremendously simplify our calculation. The results 
are discussed in section 5. In Appendix A, we will review the representation theory of S U(2 ) group, including the notation of 
intertwiners and basic components of Brink’s graphical representation and some rules of transforming graphs. The detail proof 
of some identities used in main text are presented in Appendix B and Appendix C separately. 

2 Preliminaries 

In this section, we briefly summarize the elements of LQG to establish our notations and conventions. The classical starting point 
of LQG is the Hamiltonian formalism of GR, formulated on a 3-dimensional manifold X of arbitrary topology. With Ashtekar- 
Barbero variables [24, 25], GR can be cast as a dynamical theory of connection with S 1/(2) gauge group. We denote spatial 
indices by a,b,c,... and internal indices by = 1,2,3. The phase space consists of canonical pairs (A' a ,E“) of fields on 

X, where A' a is a connection 1-form which takes values in the Lie algebra su( 2), and is a vector density of weight 1. The 
densitized triad Ef is related to the cotriad e' a by E" = | sgn(det(e^)), wherein e ahr is the Levi-Civita tensor density of 

weight 1, and sgn(det(e^)) denotes the sign of det(e^). The 3-metric on X is expressed in terms of cotriads through q a b = e' a e J b 6ij. 
The only nontrivial Poisson bracket reads 

\K(x), E b (y )} - K/36 b a 6)6\x,y ), (2.1) 

where k = 87 :G, and /j is the Barbero-Immirzi parameter. The behavior of connection under finite gauge transformations is 

A „ A s = _ (dg )g- 1 + g Ag ~ l . ( 2 . 2 ) 

The fundamental variables in LQG are the holonomy of the connection along a curve and the flux of densitized triad through a 
2-surface. Given an edge e : [0,1] —» X, the holonomy h e (A) of connection A' a along the edge e is 

h e (A) = r Pcxp (f/h4l dti dr 2 ■ ■ ■ d t„ A(e(ti)) ■ ■ ■ A(e(t „)), (2.3) 

wherein A(e(t)) = e a (t)A' a (e(t))Ti, with e a (t ) being the tangent vector of e, and r,- := -/cr,/2 (with cr ( - being the Pauli matrices), 
*P denotes the path ordering which orders the smallest path parameter to the left. Define a combination o of two edges e \, e 2 
satisfying ei(l) = e 2 ( 0 ) as 


/ ei (21), te[0,i] 

[ei ei]() ' \e 2 (2t-l), te[±, 1] ’ 

(2.4) 

and the inversion of an edge as 


e~\t) := e(l-2t). 

(2.5) 

The holonomy (2.3) has two key properties 


^eiO£ 2 C^) — he\ (,^)he 2 — h e (A) 

(2.6) 


1 A similar scheme is introduced independently at almost the same time by E. Alesci el. in [21], 
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The transformation behavior (2.2) of connection A under gauge transformation leads to the corresponding transformation behav¬ 
ior of holonomy as 

he(A g ) = g(b(e))h e (A)g(f(e)r 1 , (2.7) 


where b(e),f(e) denote the beginning and final points of e, respectively. 

Consider a finite piecewise analytic graph y in 2, which consists of analytic edges e incident at vertices v. We insert a pseudo¬ 
vertex v into each edge e and split e into two segments s e and l e such that e = s e o /~ ] and the orientations of s e and l e are all 
outgoing from the two endpoints of e. We call the new graph the standard graph obtained from the original graph by splitting 
edges and adding pseudo-vertices. Denote the standard graph by y, the set of its edges by E{y), and the set of vertices, containing 
the true vertices v and pseudo vertices v, by V(y). Our following discussion is based on the standard graphs. 

To construct quantum kinematics, one has to extend the configuration space 3K of smooth connections to the space iR of 
distributional connections. A function / on ^ is said to be cylindrical with respect to a graph y if and only if it can be written 
as / = f y o p y , wherein p y (A) = ( h ei (A ),.., h en (A)) and e\,..,e„ are the edges of y. Here h e (A) is the holonomy along e evaluated 
at A e R and f y is a complex-valued function on S U( 2)". Since a function cylindrical with respect to a graph y is automatically 
cylindrical with respect to any graph bigger than y, a cylindrical function is actually given by a whole equivalence class of 
functions f 7 . We will henceforth not distinguish the functions in one equivalence class. The set of cylindrical functions is 
denoted by Cyl(JH). The space Cyl(JH) can be completed as a Hilbert space, the kinematical Hilbert space V7 kH1 . 

Given a standard graph y, we assign each edge e with an unitary irreducible representation with spin j e , two edges incoming 
to a pseudo-vertex v with the same representation, each true vertex v with an intertwiner i v , and pseudo vertex v with the conjugate 
intertwiner it. Then the spin network state reads 2 

T rJ.M) ■= ® h ■ (g) XjShe(A)) • (g) il , (2.8) 

veV(y) eeE(y) veV(y) 


where ■ stands for contracting the upper (or former) indices of representation matrices 7tj e (h e (A)) with indices of intertwines i v at 
true vertices v, the lower (or later) indices of 7ij e {h e {A)) with indices of conjugate intertwiners it at pseudo vertices v. Given n 
edges with n spins j i, • • • , j n incident at a true v, matrix elements of the intertwiner i v associated to v take the complex conjugate 
of (generalized) Clebsch-Gordan coefficients (CGCs) up to a constant factor (see Appendix A for detailed explanation), i.e., 


(i**) 

\ Jn 


M 


= (i ] - s .) 

V Jv-Jn) r 


M 


:= (-1 ) 71 Z ' =2 - /i J (JM;d\jim l j 2 m 2 --■ j n m n ) 


= ZU J i J ^ (a 2 k 2 \j l mij 2 m 2 )(a 3 k 2 \a 2 k 2 j 3 m 3 )---(JM\a„. l k n . l j n m n ), (2.9) 


k 2 ,-,k ,,-1 


where (7/17; a \ j\m\j 2 m 2 ■ ■ ■ j n m n ) is the complex conjugate of generalized CGCs, which describe coupling n angular momenta 
j i, • ■ ■ ,j n to a total angular momentum J in the standard coupling scheme (that is, we first couple j\ to j 2 to give a resultant a 2 , 
and then couple a 2 to j 3 to give a), and so on), and a = {a 2 , ■ ■ ■ ,a„ 1 ) denotes the set of the angular momenta appeared in the 
intermediate coupling. Notice that the intertwiner presented in Eq. (2.9), differing the factor (-1 y i_ 2 i= 2 tW f rom CGCs, is more 
convenient to be simply represented in graphical formula. Matrix elements of the conjugate intertwiner it associated to pseudo 
vertex v at which two edges with the same spin j incoming are given by 

(iff'" 2 = (iT) 0 n,n2 := O'my'nilOO). (2.10) 

The assignment of intertwiners to the true vertices and conjugate intertwiners to the pseudo vertices is compatible with the 
transformation behavior (2.7) of holonomy. The CGCs are usually chosen to be real so that 


(j\m\j 2 m 2 \JM) = (JM\jimij 2 m 2 ). 


( 2 . 11 ) 


It is, therefore, not necessary to sedulously distinguish intertwiner from its conjugate when we do calculation. 

The gauge-invariant spin network states, which keep invariant under gauge transformations, correspond to the states whose 
intertwiners in (2.9) associated to true vertices are specially chosen such that the resulting angular momenta 7 = 0. The spin 
network states which are variant under gauge transformations are called gauge-variant spin network states. The normalized 
gauge-invariant/variant states consist of the orthonormal basis of the gauge-invariant/variant Hilbert space [26], 

Two elementary operators in LQG are the holonomy and the flux operators. The holonomy operator acts as a multiplication 
operator, while flux, essentially the self-adjoint right-invariant operator J 1 , acts as a derivative operator. 


[h e ,(A)] B c ■ fyihei (A), ■ ■ ■ , KSA)) := \n tl2 (h e ,(A))] B c f y (h ei (A), • • ■ , hJA )), 


7‘, • f y (h ei (A), • • • , h e ,(A), • • ■ , h en (A)) := -i - 


/=o 


f y (h ei (A), • • • , e tTi h ei (A), ■ • ■ ,h e JA)). 


( 2 . 12 ) 

(2.13) 


2 The spin network state can be normalized by multiplying \/2 j e + 1 for all e e E(y). 
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3 Graphical method for LQG 

3.1 Algebraic formula 

In LQG, under different physical considerations, one needs to construct operators, e.g., the geometric operators and the Hamil¬ 
tonian operators, corresponding to their classical quantities based on the two elementary operators h e {A) and J' e . The action of 
those operators on a given spin network state will involve the actions of the two elementary operators. The action of h e (A) on the 
spin network states involves essentially the decomposition of the tensor product representation of S U (2), which is well known 
as the Clebsch-Gordan series 


nx injm m \ 2 = X A iJ iA 


n 


(3.1) 


where (i J . ) 

V 7i 72/, 


J,M,N 

~ ( 1 y i 7*2 J{ T ATI i. 7i. 7*— 7i — \ Tlia font tbot tbo nnnrotnr jJ 

/2/ni«2 

results in 

= » 2 > 

Given a spin network state T y /fA) on a graph y, we consider a true vertex v e V(y) at which n edges e\, ■ ■ ■ e„ incident and 

he,m B C ' 


(—1 ) y 72 J (JN\jinij 2 n 2 ). The fact that the operator i J . . is unitary and its matrix elements take real numbers 

«1«2 7l72 

M 


y» 7 .* 


r v 

r r,/,L 


denote P’ ,/A) the terms, in T •■(A) , directly associated to v. Then the action of the holonomy operator [h ei (A)] B c on T v ,,(A) 


reads 

[Ji ei m B c -T" AA) = U’ a ) 

^ y, i,i \ /; 


y» 7 »* 


M 


[^(VXT SI ■ ■ • K^,)] ’ ' • ["Mr nn 


m i —mj--m n 


/ -*\ M \\ m/7) / \ 

..«. E -w.- w> 


Jr” 1 ,’"! 


where [iTj I {h ei )\ m, ni = [7r J; (/! e ,(A))] m ' ;i; . On the other hand, J' e is a derivate operator acting on a given spin network state. The 
action of J l e/ defined in (2.13) on T v ,,(A) reads 


yj,' 


.4 • r ,ta) = ( 

' y.r-i 

V ;S ) 


_ / 

’:Jj\ 

- 1 

l v J 


(i J ’ d ) 

"1 

V v / my-m'j—m, 


M 




-l — 

d t 


t =o 


[njXe n ‘he I )] mi ni )---[n j SheJ] m ’- nn 


= ^ M mi ni • • • (-i K (T,-)r ffl ;K k;<A,)l 




(3.4) 


which indicates that J 1 leaves y and j invariant, but does change the intertwiner associated to v by contracting matrix elements 
of the z-th r with the intertwiner in the following way. 


= (i?*) , (-'K(r,)]”iJ. 

1 \ / ni\ —mj—m n \ \ J 


(3.5) 


However, in practical calculation, it is not convenient to directly compute the contraction of matrix elements of r, with an 
intertwiner. One usually introduces the irreducible tensor operators [27], or the spherical tensors of r,, to replace the original r, 
(for a reason that will become clear in a moment). The spherical tensors r ; , (/j = 0, ±1), corresponding to r, (i = 1,2,3), are 
defined by 


fo := r 3 , r ± i :=+-^= (n ± ir 2 ) . 

V2 


(3.6) 


Then the contraction of matrix elements of r, with an intertwiner is transformed to that of their tensor operators with the inter¬ 
twiner. The matrix elements [7r ; (r ;; )]"' can be related to the 3 /-symbols (or CGCs) by (see Appendix B.l for proof) 


Mr,)]™' -\j2j(2j + l)(2j + 2) P j,, 


| y^m l. 

ill I D 


(3.7) 


where C“ m ' 


:= (—1 ) /+m is the contravariant “metric” tensor on the irreducible representation space 77, of S U( 2) with 

spin j (see Appendix A.l for a detailed explaination for the C"'” m ) [28], The spherical tensor r ;j generates the self-adjoint 
right-invariant operator 7^ defined by 


Je, ■ fy(h ei (A), • • • ,/z C; (A), • • • ,/z e „(A)) := -i 


dr 


f=0 


fy(h et (A), ■ ■ ■ ,e‘^h e ,(A),--- ,/z c „(A)). 


(3.8) 
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The action of ./£, on 


(‘n 


M 


reads 




J e, 


{?*) * = (?*) , 


(3.9) 


Any gauge-invariant operator, e.g., the volume operator considered in this paper, defined by J's can be expressed in terms of the 
corresponding J M s. Hence its action on the spin-network states essentially is equivalent to contracting 3 /-symbols (or CGCs) 
with corresponding intertwiners. 

3.2 Graphical representation and graphical calculation 

The basic components of the graphical representation and the simple rules of transforming graphs are presented in Appendix 
A.l. In graphical representation, the 3 /-symbol is represented by an oriented node with three lines, each of which represents a 
value of j, namely 


J l J 2 J 3 
m\ m 2 m 2 


(3.10) 


where - and + denote a clockwise orientation and an anti-clockwise orientation respectively. A rotation of the diagram does not 
change the cyclic order of lines, and the angles between two lines as well as their lengths at a node have no significance. The 
“metric” tensor C^? in Eq. (A. 7) which occurs in the contraction of two 3 7 -symbols with the same values of j is denoted by a 
line with an arrow on it, i.e.. 


c 0) =(—i v~ m s , = (~i y +m 's ’ = 

^ u m,—m V ± J u m,—m m - 

and its inverse in Eq. (A. 8 ) can be expressed as 

c mm ' = (_i y~ m s , = 

'-'(j) V 1 ) u m,-nr V u m,-m' m - 


(3.11) 


(3.12) 


Summation over the magnetic quantum numbers m is graphically represented by joining the free ends of the corresponding lines. 
The contraction of a 3 /-symbol with a “metric” is represented by a node with one arrow, which provides us a way to representing 
the CGC, e.g.. 


(j3m\jimj2m 2 ) = (- 1) J1 h 13 yjlh + 1 

= (-i y'- h ~ h V2 h +1 


(■n 

ji 

J3 ' 

j C (7s) 

\m i 

m 2 

m 3, 

mi 

Jl 


■m 3 = (—1 y" 1 








m 2 





1 - 72-73 


mi 

•M, +1 '■[_ 


V / n 


(3.13) 


in Eq. 


This is the main motivation for Brink to modify the original Yutsis scheme [29, 20]. Hence the intertwiner ( 

(2.9) associated to a true vertex v, from which n edges are outgoing, is represented in graphical formula by Eq. (A.47) as (see 
Appendix A.2 for a detailed interpretation) 


(*v /;a ) M = {ij’ a j ) M = h V 2 Oi + 1 V27+1 ■ J 

V / m\m 2 "-m n Wi Jn/m\m, 2 "-m n 11 V 

i=2 


(3.14) 


Now we will extend Brink’s representation and propose a graphical representation for the unitary irreducible representation TXj of 
S U(2). The matrix element [nj(g)] m is denoted by a blue line with a hollow arrow (triangle) in it as 


[nj(g)r n =■■ m. 


!>■ 


(3.15) 


The orientation of the arrow is from its row index m to its column index n. The Clebsch-Gordan series in (3.1) can be represented 
by 


mi - 


J i 


m 2 - 


32 


A ^ 

= Y^h + 1 ) 


mi 


J 1 +J 3 


P 2 

m 2 


J3 


J3-J1 


J1+J3 


72 1 

m'n 


< = J](2j3+ l) m '^ 


r 2 

m 2 


3 3 

4>- 


J3-31 


J2\ 


(3.16) 
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where, in the second step, we used (A.52) to flip the orientations of two arrows. The matrix element [nj(g )]" can be presented 
in Eq. (A.66) by 


[xte- 1 )r m = cZ'[n j (g)] m ' n ,Cti= n - 


Using Eqs. (3.17) and (3.16), we have 


7 i 

m 1 -p£>- 

h 


■ m 1 m i ■ 


m 2 - 


■ m 2 m 2 - 


31 

S- 

h 

4V 


& 


4 


(3.17) 


3 3 

Z m-. •^ 1 +-Z 3 hSv 33—3l 

< 2 A + 1 > - p *f 


v-, 71+73 IN. 73-71 , 

p —&>— -p-. 

m 2 


73 


73 


V* 

ro 2 


(3.18) 


For group elements as the holonomies h e = /z,,M) of connection A along an edge e, their matrix elements can be simply repre¬ 
sented by 


Me)]\ = 




(3.19) 


By Eqs. (2.6), (3.17) and (3.19), the matrix elements of the inverse of a holonomy can be represented by 


[*j<h; 1 )r m = = «- 






-<3- 


-771 = n ~ 


3 

Eqs. (3.16) and (3.18) yield the coupling rules of representations of holonomies as 

e / 


7711 - 
777 2 - 


3 1 
e 


7771 


J 2 


and 


7771 - 


Jl 

^-1 


7771 


777 2 - 


32 


Then the spin network states T v ^(A) associated to v can be represented by 

y,j,i 

r AA) = (i v J ’ a ) M [n h (hjr ... [:Kj,(h ei )] • ■ • [njMr- n 

y,j,i \ v tm\—mi—m n J «i Jl 1 n i Jn n n„ 


jq V 2 ai +1V27 +1 

7=2 


■ M 


«/-i a/ 

The action of [/z e/ (A)] 5 c on the spin network state given in (3.3) can be represented by 

nj C 


n -1 


\h ei (A)] B c ■ r y ;f A ) = fl V2a,- + 1 V27 + 1 

i =2 


'A 

‘ji e i‘ 

l 

1 2 e ni 

> jn 

jl 

V ^ 

ji 

7 

jn 


a-i 

a n _i 

J 


■ M 


(3.20) 


71+73 

e 

33-31 

■ mi v—1 mi 

- 2(273 + 1 ) 

7 l +73 

e 

73-71 , 


\j2 

73 

u 

|j 2 

73 

-►“!- m i 

7 2 | 

( 3 . 21 ) 

777 2 


m' 2 

73 

777 2 


m' 2 


ji+h 

e 

33—jl 

mi 1 ™.i ■ 

= 2(273 + 1 ) 

71+73 

e 

33—jl 1 


fa 

33 


k 

73 

- ►“!- m 1 

( 3 . 22 ) 

777 2 


m' 2 

73 

777 2 


m 2 



(3.23) 
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n— 1 

= Y] V 2 Oi + 1 V 2/+1 £( 2 /, + 1 ) 

*= 2 J'i 


n, 


3i 



— 

^— (' 


A 

k 2 



‘A 

l n 

,ji i 

, 1 

2 

-*4-B 

L 3 n 

3i 



3i 

3i 

jn 

\ ...» 




a i-i cLi a n -i J 


(3.24) 


The spherical tensor in Eq. (3.7) can be represented graphically by 


[nj(= l -z V2;(2 ./+l)(2; + 2) 1 = ^ ^2j(2j+Y)(2j + 2) 

A m ! ——:— m A 


(3.25) 


U ( J' (1\ M 

Hence the action of K on i.,’ ) in (3.9) can be represented by 

1 \ J m^m-y-mn 


A ■ = \ V2;(2; + 1)(2; + 2) f| V^T VtTTT 


n-1 


i=2 


X e. 



Jn 


. (3.26) 


■ M 


The graphical calculation is to use the simple rules of transforming graphs (see A. 2) to derive the result. 

The starting point of our scheme is the so-called standard graph y st d, which is obtained from its original graph y org by splitting 
edges and adding pseudo-vertices. We still need to show that the spin network function associated to the original graph y org is 
equivalent to the one associated to its corresponding standard graph y stc i when acted by an operator, e.g., the volume operator. 
Recall that the standard graph y st d is obtained from y org by the following procedure. We insert a pseudo-vertex v into each edge 
e of y org and split e into two segments s e and l e , such that e = s r , o ' and the orientations of ,v ( , and /,, are all outgoing from 
the two endpoints of e. The standard graph y st d consists of the new segments .v,, and l e , the new adding pseudo-vertices v, and 
the (old) vertices of y org . We can transform the spin networks based on the original graph into those on its standard graph by 
explicit transformation rules, and then find out their relation. Consider an edge e with irrcpresentaion j in y org starting from v 
and ending at V , assigning the intertwiners i v and i V ', respectively. We assume the edge e in the original graph y org is regarded as 
the k-th edge and the k'-th edge in the set of edges which incident at v and v' respectively, i.e., b(e) = b(ek),f(e) = f(e' k ,). The 
relevant ingredient of a spin network state associated to the edge e takes the form (see (A.48) for the graphical representation of 
the intertwiner associated to v at which there are coming and outgoing edges, (3.19) and (3.20) for the graphcal representation of 
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the holonomy) 




,/ -/ - m'r -M' 

ai'*).- ‘ 




H^ 1 x/2«|+lx/27 + 1 ■ 


Mhe] 


(3.27) 


a fc-l <2 fc 


The transformation from y org to y std induces the following transformation in the spin network state. 



(3.28) 


where we have used (3.20) in second step and the rule (A.50) to remove two arrows in the last step. Repeating the above 
procedure, we can transform the spin network states associated to the origin graph into those corresponding to its standard graph. 
By this trick, we finally find out the corresponding relation of the spin network states between the original graph y org and its 
standard graph y std . The intertwiners associated to the vertices of the origin graph is replaced by its standard formula, and the 
adding pseudo vertex v for each edge e with a di-valent intertwiners is graphically represented by an arrow with orientation 
opposed to that of the original edge 3 . 

Let Y' vei be an operator assigned to vertex v and edge c/ intersecting v by the following formula. 


f A’.e,;/. ’ fyJK > ' 

> hei) 

:=i /roig^,," 

• ,he,e ,T < 


when ei is incoming 

|/v,£',;k ’ fy^ ei , ■ 

5 h ei ) 

• = $l f =0^org(^l.- 

•• ,e tT ‘h e „- 

? hej\ 

when e/ is outgoing 


(3.29) 


Our task is to show that the action of L on an arbitrary given spin network state of y org is equivalent to the action of J\, t . R 
on the spin network state of its corresponding standard graph y std . The proof can be done by using their tensor operators J 11 
(jj = 0, +1) in the graphical calculus. Notice that 


^y^Ahe)r m . = j%, L -[n j (h e )r m . = ij t 


t =0 


[nj(h e e ,T ")] m m , = i - 


t=0 


{[nj(h e )] m MAe'^r m .) = WhaTnUnfa)]” 


(3.30) 


3 The intertwiner, a line with an arrow, associated to the pseudo vertex v is not normalizable since we adopt nj e (h e (A )) rather that its normalized form 
1 7ij e {h e (A )) in the spin network function (see Eq. (2.8)). If the original spin network is normalized, the intertwiner associated to v will automatically be 
normalized. Then it will be expressed as times a line with a spin -j arrow in the graphical representation in Eq. (3.28). 
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We obtain 


'Xir=2 1 \/2oJ + iVaj' +1- 


f- 

v',e\L 


11 2 "t 1 \/‘2J + 1 


o! k , a'k'-i 


a k _i a fe 


'fli=2 >/ 2a i “I" iv^FTt- 


<v - a ‘-' 

—4— 


| N /2i(2j + l)(2j + 2) 


flSU 1 \/2a 7 ; + W2J + 1 ■ 


jfo) 


flfc-i a fc 


■|Ti=2 + 1-• • 


a., a,./ 

- 


A^ 2 j( 2 j + l)( 2 j + 2 ) M- 


rE v‘- a i +1^57+1 ■ 


Jfe' = J 

J 




l^[ i _2 \/2a£ + lx/2.7' + 1- • • 


- ^ 


v',L; 


where, in the second step, we have used 


rJfc' _ J 

j 


itj y/2a,i + W2J + 1 • 


a k 

) 



a' k , 

a k'-i 


•jk' ~ j 

1 

j 


'j 

le 

'j 

S e 

‘3 


jk = j 

a k - 1 

a k 




a ., _ <V 

-*#- 


a*-! a fc 


Jfc' =J 

J 

j 


dk-i a k 


(3.31) 


(3.32) 


here we have used (3.20) in the first step, and (A. 37), (A.50) and (A.52) in the second step. By the above procedure, we can also 
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transform other edges into their corresponding standard edges. Thus we obtain 


(j — (I 

i \ 




(3.33) 

y org ' ■ /std 

The above calculus and analysis complete our proof. 

4 The matrix elements of volume operator 

There are two versions of volume operator in the literature. We only consider the volume operator defined in [10, 11], which 
passed the consistency check in quantum kinematical framework and was used to define a Hamiltonian constraint operator in 
LQG [30, 31, 32]. The volume operator acts on a spin network state as 


y. T 


,j/A) E A 8 x 4 2 

veV(y) \ I<J<K,e I ne J r\e K =v 


?(<?/, e } , e K ) q IJK 


■ T 

yj.i v 


(4.1) 


where g(ei,ej,e K ) = sgn(det(e/(0),ey(0),e/s:(0))), and 


qiJK ■- -4 i£ijkJ , e I J 1 ejJ k e K ~ 


4 \Wi,SikJ l ej J k eK \ ={166uASMri - 16 Wifaf'S") =: - tf£ K> ) . 

(4.2) 


Here J 1 defined in (2.13) is the self-adjoint operator of the right-invariant vector field on the copy of S 17(2) corresponding to the 
7-th edge, satisfying 


[7 g; , J ej ] — i6ijkJ ei 6 eit ej . 


(4.3) 


and 5ij := -2tr(r,T / ) is the Cartan-Killing metric on SU(2). The action of volume operator is local, in the sense that its action 
is a sum on independent vertices. Therefore, we can focus on its action on a single vertex. The fact that the pseudo-vertices as 
endpoints of edges are divalent and the self-adjoint operators J 1 act only at the beginning points of ej implies that the summation 
in Eq. (4.1) is only over the true vertices v of y. 

Eq. (3.4) reveals the fact that the operators quK and thus V only change the intertwiners i when evaluated on T^/^A). The 
operators q IJK acts on an intertwiner by contracting the corresponding matrix elements of r, with the intertwiner. Note that 

SijJiji, ■ (i v jj ) M = (i^ d ) , , M (-K(T ; )] m ; m 8 ij [n h {Tj)T' J J 

J c/ J \ / mi—mi—mj—m n \ / mi" m I —m J --m n \ J ‘ JJ J / 


= , , [nj^T^ m, C m bry,(T v )]" 


~ _c (n f r ■ (i J ’ d ) 


M 


where in the second step we have used the following identity (see Appendix B.l for proof) 

Of' m 6 ij [KjjiT j)] m J mj = [nMV'rnC^ [n jj {T v )] m ' J 
Hence the operator q^j^' ,J> and q^j^’ JK> in (4.2) can be represented in terms of by 

ii J v J ) M = 16 C^J p ej r eK C^r ei J v ej ■ 

\ v / m\— mj—mj--m K --m n P a J eK P v 1 ej \ x / m\ — mj—mj—m*- 

= 16(i/’ ) , , , l /i jjv i vj\ ■ l /i ik 


(4.4) 


(4.5) 


Qijk 


M 


rti\ —mi—ntj—mK—ntn 

[^(r v )]<„,,Cft [Xj K (T V ’)] m ' K x K(r,)] m ; ffl K- y (v)] m G 


(4.6) 


qf// K> ■ (i dJ ) M = 16 


- 16 (i/ ; *) 


M 


m' 




[xMr- C^ [nj,(Tu')] J x [7r,(T v )]^,„„C , [ 7 tUtv'W 


(4.7) 


/ j. M 

With above preparations, we now turn to the action of q IJK on the intertwiner (L ’ associated to a true vertex 

v in the graphical method. We first consider the case I > 2 and K < n, where cq \ and a K will appear in the final result. The other 
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special cases will be dealt with later. According to Eq. (4.6), the first term in the parenthesis of Eq. (4.2) evaluated on intertwiner 
(3.14) can be represented by the following graphical formula (we present only the parts of the graph of the intertwiner which 
closely connect to the key steps in the following calculations), 




n— 1 

= q^K IJ> ■ fl V 2 «/ + 1 V2/+ 1 


nij m J+ 1 


ttlk-i m K 


i =2 


_ a,i _ a/ + i • • • a,j. 


- aj+i ■ ■ -clk-i — Ok-1 — 

m,j rrij+i rriK-i 


n— 1 


- X(jl> jj) 2 X(jj 7 jx) 2 Y\ V 2 ^/ + 1 VZ/TT 


i=2 





jj 





ji 


jl+l 

jj-i 7 

7 


jj+i 

Ik -i 

u 




r 

+ 



r 

+ 

ji' 

\ 


\ 

'jj 



f 

'jK 


— id/ _ a/ +1 • ■ • aij_2 — dj. 


aj — dj+i • • • d/r-2 — a*r-i| — 


(4.8) 


where X(j\, 72 ) = 2 7 * 1(271 + 1)(2 7 *! + 2)2 72(27*2 + 1)(2 72 + 2). Similarly, the second term in the parenthesis of Eq. (4.2) acting on 
the intertwiner can be expressed as 


"UK 


■ (■-'') 


M 


mi-mi—mj—mg—m n 


n— 1 

= X(jf,jj)2X(jj,jf{)2 J - [ y/2cij + 1 V 27 + 1 

i=2 



(4.9) 


It is obvious that the two terms and q^ JK> in (4.2) are gauge invariant. Hence the operator q IJK and the volume 

operator (4.1) are gauge invariant. Each of them leaves the intertwiner space TfY . with the intertwiners as its orhtonormal 
basis, determined by the given j\, ■ ■ ■ , j„ and the resulting angular momentum J, at the vertex v invariant. Therefore the action 
of quK on an intertwiner can be expanded linearly in terms of intertwiners in 7T‘ | . at the vertex v. In graphical language, they 

leave the vertical lines denoted by /, and the last horizontal line denoted by J invariant, but change the intermediate couplings 
ai labelling the intermediate horizontal lines. Hence, in graphical calculation, our task is to drag the endpoints of the two curves 
with spin 1 down to attach the horizontal lines, and then yank them away from the horizontal lines following the simple and 
rigorous rules of transforming graphs presented in Appendix A.2. Certainly, there are many alternative ways, corresponding 
to different choices of recoupling, to remove the two curves with spin 1. The results obtained from different ways are related 
by unitary transformations. In what follows, we choose a way guided by the simplicity principle that the number of changed 
intermediate values a, is as little as possible and the final result is as simple as possible. The calculations of the action of q IJK on 
the given intertwiner in the graphical method consist of the following four steps. 

The first step involves dragging the two endpoints of curves with spin 1 attached to lines with spins jj and jx respectively 
down to join with two horizontal lines denoted by spins a/ and u K i. To do this, we use the following recoupling identities (see 
Appendix B.2 for proof), 





' ji 

ai -1 

a 7 

1 

3k 

+ 



' 3k 

Q-K-l 

aK 


fa i- 1 

ji 

a 7 


Ji 

l 1 

a'j 

ji\ 

a/_i 

a'i 


Yp b K-i + l)(-l) aK ~ b x-' +jK+1 


fa K 

3k 

ax-\ 



Jk 

u 

b K -1 

Jk 

aK- 1 

b'x-i- 

1 

a K 


Then we have 

~<JK,IJ> 

"l J K Vv 

= Yj- 2a 'i +1 )(-i/" 


a i -1 

i 


ji 

a', 


YP- b 'K-x + 1 )(-D a 


-b'^+jK+l l a K JK 


1 K 


K -1 


a K -1 

h 


(4.10) 


(4.11) 
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n— 1 

xxUjjjyxUj.jK) 2 J - [ V2a,- + i V2y + l 

i=2 



(4.12) 


Let us turn to the second step. We first consider the case that J > I + 1 and K > J + 1. The other cases will be handled later. 
We move the two points labelled by ( a ! 1 , fl/, 1) and (clk-\ , b ' K] , 1), step by step, to the right hand side of (ay_2, jj-i,aj-\) and the 
left hand side of (a j, j j+\,aj + \) respectively, by repeatedly applying of the following identies (see Appendix B.2 for proof). 


Ypa'i + a ' a i 


— = YP b 'm-\ + l)(-l)^ +a "- 1+1 (-l)^-i-;’ 


1 Cl\ fl/-lj “l-i " “i 

Jm b' m -\ K 


(4.13) 


(4.14) 


where l = I + 1,- ■ ■ ,J - l and m — K - 1 ,■■•,/ + 1. Then we have 

Y Y ( 2 «/+i)(-i ) a,_i ' a ' +j/ ^ a/ “ 1 jl ai 

( a r a i+i'" ’ a j-0 Wk-v'" >bj+i>bj) 

J -1 


r 2 , };}x(2^ 1 + D(-ir 


b' K _i+jx+l l a K JK O.K -1 


1 6' 


■/* 


x n ( 2 a; + i)(-i) a ;-i +a '-' +i (-i) a '-'- fl ' + ^ J^ a i -1 

II I 1 a/ a/_i 


/=/+i 

a:-i 


n + Df-D 1 


Z/ 1 +6hjj—1+1 


(-ir 


m—J+l 


_,~jm ba 

1 Cl m d m — 1 


, W 

m -1 m 


nil mi+i 


TTij—i m j 


rriK-i m K 


xXUiJjY-XUjJk) 2 ]~[ yj2cij + 1 V27 + 1 

i=2 


jl 

ii+i 

jj 

jj -1 

jj 

A 


jj+i 

Ik- 1 



f 

+ 

b', 


- L 


- a /+1 •• • aj_ 2 ~ ~I 


gXJ 


(4.15) 


where ( a' I ,a' I+l ■ ■ ■ ,a' J X ) and (b' K V ■ ■ ■ ,b' j+l ,b'j) denote the all allowed tuples, and aj and b' m l are determined by a' l _ x and b' m 
respectively. 

In the thrid step, we dragg the two endpoints of two curves with spin 1 attached to the line with spin / / down to join with two 
horizontal lines denoted by spins aj-\ and cij, respectively. To do this, we use again the identity (4.11) to pull down the point 
( jj,jj , 1) and simplify the two quantities nfj/ + i( _ iy , '~ ,_ “' +J ' and Y\m=\+ \ (~ I _J ' n . Then we obtain 

Y Y (2a'/ + l)(-ir-'- a '' +J '\Y' j a/ |x(2^_, + lX-lfr-^-i+A+i/ a * J* Uk ~' 


l ' a 'r“’u i 1 ) 0>'k~ p" •bj+iJ’j) 


1 b' 


JK 


J -1 


FI (2a! + l)(-l)°;-i +a '- 1+1 P' a '-' a i 

1 1 II d[ d[-\ 


1=1 +1 
K -1 


m=J +1 

yij-b’j^+jj+l J a J J] ay_i 


; L Ff + l)(-l)*L- 1 +«™-l+ 1 {■/» ^m-1 

,1, (1 a;;; d,. : 

Ypb'j-\ + 1)( — 1) 


1 1 


my my+i 


m K -\ m K 


n -1 


x ^ 07 > yV ) 2 yV ) 2 ]~[ ^2dj +1V27 + 1 

i=2 


h 

ii+i 

jj 

jj-i 

jj 

A 


Jj + 1 

Ik- 1 

JK 



aj~i b'j j 


fe'y 







. 


(4.16) 


- &:, + i ■ • • b' K 
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The identity in Eq. (4.10) can be used to drag the remained endpoint, which attaches to the jj line, down to join the horizontal 
line denoted by aj and yield 

Y Y ( 2 «'/+ 1)(-1 )“'-■-W a '-’ J j, a/ | X (2b' K _ x + l)(-l) asr_ ^-i + - /sr+1 [ Uk A k Uk - 1 




1 b K , Jk 


y-i 


Ff (2a; + l)(-l) a; - 1+a '~ 1+1 H' a, - x a1 

t 1 1 I 1 ci i ai -1 


/=/+1 
K -1 


X (_l)a/-ay-i+Z l=/+ i7i 

x(-l)^-i-^-2-«^ Ff (2^' )I _ 1 + l)(-l)*™-i +a ”- 1+1 1 7 !" /?m - 1 

, , II 1 

m=J +1 v 

x Y (2b 'j-i + 1)(-1 +J>1 |7 b i J a (“'l x ^](2aj + 1 )(-]/)-! J °>-i 2y 


b'j_\ jj ay 
y 2 'y 


x X(ji, jj) 2 X(jj, j K ) 2 J””[ ^20/ + 1 V27 + 1 

;=2 



jr 

7'y+i 

7-y-i 

a/,! 6j, 

jj 

a'j tfj 

jj+i 

7k-i 

7 k 

* “. ‘-trcr .^. ■ 


■ (4.17) 


Let us go to the forth (last) step. It is the time to remove the two curves with spin 1. To do this, we use the identity (see 
Appendix B.3 for a proof) 


,, a 'j ~ aj ~ u 1 , On! W 

- 1 4 ) * 3 = (-\) a l~ aj+X - /’/, fc', 


2 aj + 1 


(4.18) 


to remove the two curves with spin 1 from intertwiners and obtain 


Y Y ( 2 A/ + i)(-i) a '- | - a;+J '<r / “ 1 J j, a ' i x (2b ' K _, + !)(-! f K - b '«-' +jK+i -j aK ~ l 


7-1 


_i \fl/-ay_i+Xt/ + i jl 


X (-1 )■ 


n (2«; + ix-ifi-'*"* 1 f ^ 


/=/+l 
AT—1 


X (-i)^-.-^-2S + iy m rf (2^_ 1+ 1)(-1)^.— +1 fr ^ 1 

i 1 II d m d m -1 I 

m=7+l v 7 


: ^(^-t + 1 )(-D“ 


j_ 1 +y'j+i J a i Jj fly-i 

1 b'j_ 1 jj 


Yj2aj +1 )(-i ) fe (-i _a / + T ' J j y 


x(-l)' 




2(7' , + l^->’^ X( 1)1 




*7-1 


2 aj + 1 ' 


X 2 f(j/, jj) 2 X(jj, j K ) 2 |~[ V2fl; + 1 V 27 + 1 


m/ rm+i mj-i mj "‘j+i m.K -1 

7/ 7/+i Jj-i jj jj+i jie-i 7 k 


7=2 


- a', - Oj +1 • • • Oj_ 2 - aS_! - 6j -b' J+1 ---b' K 

Summing over b'_. and aj and relabeling the indices b' by a', we get 


(4.19) 


E (2aj + 1)(-1 j 7 . -1 J \ a1 

(a' r -,a' K _ j) 1 1 a i JI 

-a' K _ l+ j K+ \ fflA- 7 a flA-1 


x( 2 aj. ! + 1 )(— 1 ') clK ~ a K-\ + Jk + * ) 

A 1 I ^ ^Kj-l •/£” 


from the first step 


X (-i) a '-^- i+2 w+i-'' n (2a; + i)(-i) a ;-i +a '-i +i -j 7 / a, -> a ‘ 


/=/+1 
a:-i 


1 < 2 / £f/_ 


x(-l)^-i- a (-2™:) + iy™ n (2a'_, + l)(-l) a m-i +a ”- 1+1 


m=7+l 


1 dm dn 


from the second step 


13 













































aj J f j aj ~ l L (_l)a/-flj-l + J> +1 (-l)^ a J-l +fl -'-i- 1 

1 a j-\ jj J 

x i Uj -1 p aj l (-1 ) a ’j-r a ’j + ]j(-]y'j-aj+i 

1 a , i.r\ 


from the third and forth steps 


x X(ji, jj ) 2 X(jj, j K ) 2 J - [ y]2di + 1 V 2 / + 1 


m /+ i 

U i i 


m K _i m K 


Ik-1 Jk 


a, ~ aj. 


Up to now, we have completed the key calculations and obtained the action of the operator qj J j^’ IJ> on the intertwiner ( 
as (4.20), which is a linear combination of new intertwiners. 

Now we want to simplify the expression (4.20) under the following two considerations. One is to get more symmetric factors 
in the two multi-products. The other is to simplify the exponents. Notice that the result (4.20) was obtained in the case of / > 1 +1 
and K > J + 1. Under this case, the product terms of (2a' + 1) in Eq. (4.20) reduce to 




(4.20) 

M 


K-\ 


K-\ 


( 2 a' K _ { + 1 ) f[ (2a ' m _, + 1 ) - (2aj + 1 ) [ ( 2 a m + 1 ), 


(4.21) 


m—J +1 


m—J +1 


which enables us to write the multiple product over m as the formula which is more close to the multiple product over / in Eq. 
(4.20). Thus we have 


Y (2a'i + 1 )(- 1 J a '" 

m 'k I * 

J-l 

x ji 


cii -l Ji ai 
a', ji 


x (2a'j + 1)(-1 ) a 


a’ K _ l +j K +1 I a K JK Cl K -1 


1 a' 


K- 1 


JK 


r\(2a’ l + l)(-l) a ^ +a ‘- 1+l \ J ‘ a >-> a i 
ii I 1 ai ai-i 


i=i +1 
K -1 


x(-l]““[ (2a' m + 1 )(—l) 0 "- 


m=J +1 


+ 1 J Jm Cl m _ 

1 a n 


dm -1 


fly Jj fly-i 


1 a’, 


y-i 


Jj 


(- 1 )' 


aj-a'+jj+l 


(-i r 


a' ,+flj-i-l v I Ay -1 JJ Ay 


1 a', jj 






xXU^jj+XUjjKP f] V2fl; + 1 V2/ + 1 

i =2 


m, m I+1 

)i +1 


(4.22) 


The exponents in (4.22) can be simplified as 

(_n a2 -i _ ^ + y'(_l) a2 ^ i 4^ + y^ +1 ('_n a22;fl ?=r+2]£>+i yir_^y>l^>c a )-Em=/+i iY'j-i~ a 'j + J J (— ly'jxfr 1 

- c-iyi-i+ji+an+iKC-iyi-aW- i)2tAiy'i(_i)-2mw+iy™(_i) a )-i +a )(_i)- 2 (‘4 + ‘4-i-fd 

= (_i) a i-i + y/ +fl « + y«(_i)° 2 -a j(_i)2£/liy/(_i) - 2^/+iyi"(_n o y- 1 +a / ; (4.23) 

where we have used, in the second step, the fact a'j + a' JX - jj e Z due to the allowed triple (a'j.a'j { ,jj) satisfying triangular 
condition. Considering the simplified factor and properly adjusting the ordering of multi-products of V 2 a + 1 , we finally have 
the compact result 


& J> ■ {#*) 

-- Y (-i ) a, -' +ii+aK+iK (-\) ai -*+-\)^ h (-\r ^ j -x(jj,jj)*x(jj,j K p 




7-1 


ffl/-l 

7/ 

A/] 

ja K 

)K 

l 1 



\1 

a K- 1 


JK 


* n V (2 ‘ , i +i)(2 “ |+i)<-i) ' i '""*‘ii z, 


l=l+\ 
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m—J +1 


7-1 


2 m + 1)(2 a m 

+ l)(-l) fl -l +a " 

.-i+i. 

I Jm 

I I 

jaj 

j] 

«/-ll 

. K-i 

Jj 


u 

a ’j -1 

jj 1 

l i 

a) 

jj) 


K-\ 


n -1 


] V2a,- + 1 Y\ V 2 < + 1 V2a* + 1 V27 + 1 


d m -\ 


m j mj +1 


rriK-i m K 


i =2 


S—I 


k=K 


— a I — a I+1 ■ ■ ■ a,. 


- a J+l ■■ ■ a K 


(4.24) 


The above result is in the case of / > I + 1 and K >./+!, which assures the limitations of the summations and multi-product 
exist, i.e, the upper limitations are always not less than the low limitations. Intuitively the results for the other cases can be 
seen as special cases of the above result by omitting the corresponding summations and multiplications which do not exit. The 
following analysis in step by step will confirm this intuition. 


(I). The first case of / = I + 1 and K — J + 1 

In this case, the corresponding processes can be simplified and shown in the following way 


mi m,i m K mi m.i mK 


ji "\ U 

( \ ^ 

ji 

ji ~ -\ j K 

. * - 

"-■g-g-" 

g g - 

mi m.i mK 

ji Jj Jk 

aj b'j-.a’j 6' k , 

m, mj m K 

jl jj jK 

1—> 

g.g 

-1 „r - j: - 


The result can be directly obtained from (4.20) by omitting those terms appeared in the second step as 


V (2a'j + \)(-\) a, -'~ a ’i +il J a/_1 J \ a, \ x{2a’ K _ x + l)(-l) fl ^-i + ^ +1 \ a f \ K aK ~ x 

la!, a' K _,) ( 1 a i Jl ) l 1 a K -1 Jk 


a J Jj Oj-\ 
1 a'j_ 1 j] 


. (_ 1 yj-a-'j -1 +//+ 1 (_ 1 y“'j -1 +«/-1 -1 x a J 1 

1 1 a!i i.r I 

1 

x X(ji, jj) 2 X(jj, j K ) 2 J””[ y/2a, + 1 V2 J + 1 


aj +1 


n -1 


mi mj 


1i l 


i =2 


(4.25) 


(4.26) 


The exponents can be simplified as 

— (■_ ^ +ji+°K+jK(_i yi-a'j ^ y'j_ i+a'j 


(4.27) 


where, in the second step, we have made the replacements aj-\ = ai and a' K _ l = a'j. After replacing (2a' K t + 1) by (2a) + 1) and 
properly adjusting the ordering of multi-products of V2a + 1, Eq. (4.26) reduces to 


y (-iy'-' +ji+a « +j *(-ir- a ‘XUi,jj)' j x(jj,j K y- 

(■ a i" a K- 1 > 


^(2aj + l)(2a/ + 1) J(2a) + l)(2aj + 1 ) 


a i -1 
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ji 

a/1 

[<Jk 

Jk 

ClK- 

a 'i 

j/f 

u 

cl' 

U K -1 

Jk 


15 



























































J aj 

Jj 

ay-i) 

,f a 'j -1 Jj a A 

u 

a J -1 

jj j 

\ 1 «y Jj) 


/-I K -1 n -1 

x | | ^Jlcij + 11 | -yj2,d' s + 1 | | ^Icik + 1 V27 + 1 

i=2 5=/ A:=AT 


(4.28) 


As a special case, the result (4.28) can also be directly written down from Eq. (4.24). In (4.24), the multi-products FI/wVi an d 
n *=} +1 and summations Efjz+i an d XmZj+\ do not exist for J = I + \,K = J + \. Thus the terms involve these multi-products 
and summations can be omitted, and this yields the result (4.28). 


(II). The second case of / = I + 1 and K > J + 1 

The result can be directly obtained from (4.20) by omitting some terms appeared in the first line from the second step. Taking 
account of Eq. (4.21), we obtain 


X (2a, + l)(- 1 ) a '- |_a ' +/ ' j a/ ~' J j f "'j x (2aj + l)(-]/'‘-^-i + ^ +l | 


flK-a’^+jK+l \ a K Jk a K - 1 


’ a K-0 


1 dj, 


JK 


K -1 


x(-l)<-.- fl (~2,t) + iA Ff (2a;„ + l)(-l) fl - +a '"- 1+1 7 r am 

11 II d m @m- 1 

m=J +1 v 

x _)ay Jj a J-^\^_iyj-a'j_ l +ji+\^_Yya'j- i +aj-i-^ x ) Cl J-\ dJ i y'j-aj +1 

1 a'j_ i jj J l 1 a 'j jj) 


n -1 


xX(jiJj) 2 X(jjJ K ) 2 ]^[ V 2di + 1 V 2 J + 1 

i=2 


h 1 


(4.29) 


The exponents can be simplified as 

= (_i yi-i+ii+aK+iK Z*:)+ 1 im (_ ^ y'j-i +a j (_ j y 2 ^ a 'i +a 'j- 1 -ii) 

- )«)-!+«) ; (4.30) 


where, in the second step, we have replaced a y_i by aj. After properly adjusting the ordering of multi-products of V 2 a + 1 , we 
have 


2 (-l ) a/_1 + d ,+aK+ j K (- \) a ’~ a 'i (_ 1 ) _ i- X (jt ,;y)5X0 y , 7^)1 


x j(2a' + l)(2a, + 1 ) ^( 2 a) + 1 )( 2 aj + 1 ) j ^ 1 J j, j 0 / 
£-1 

x ]"~[ V(2«;+ l)(2« l „+ !)(-!)< 


a K Jk a K _i 


Jk 
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K -1 


m=J+l 


/-I 


,+a m _i + l J Jm a 

1 1 


j a J 

j] 


. K-t 

jj a A 

l 1 

a 'j-\ 

j] j 

l 1 

a 'j Jj) 


AT—1 


n-1 


[~[ V2a, + 1 P| V2fl' ? + 1 P| V2«t + 1 V27 + 1 


i =2 


S—I 


k=K 


(4.31) 


Again, the above result can also be directly written down from Eq. (4.24). In (4.24), the multi-products n/~/+i and summations 
Y/i=i+\ do not exist for J = I + 1. Hence the terms involve these multi-products and summations can be omitted, and this yields 
the result (4.31). 
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(III). The third case of / > I + 1 and K — J + 1 

The result can be directly obtained from (4.20) by omitting some terms appeared in the second line from the second step as 


V da'j + r 1 .' J i a ! \ x (2a' K _ | + ])(-] )“‘-4 -,+/^ 1 J “ K l K ° K -> 

{ 1 a i Jl ) l 1 ' lK 


"K-l 


7-1 


_1 \ ai ~ aj - l+ ^l=M Jl 


x (-1) 


n(2fl; + ix-i) a '- i+a, " i+1 r 1 / a ’-' a ' 

il II a j @i-\ 


i-i+i 


aj Jj a J-\ (^_ 1 -)fly-a)_ 1 +7y+l r _ 1r a)_ 1 +a / -i-l * \ a j-\ Jj a }\((-U a 'j~ aj+l 


(-D 


1 jj 

xX(ji,jj) 2 X(jjjKr* [ V2 fl / + 1 V2/ + 1 


1 a'j jj 


; u_if/-i-^(-i)“ 


n— 1 


i=2 


Wj m K 


lJ Ik 


(4.32) 


The exponents can be simplified as 

= (_l) fl '-l + f' +a A + fA(_l) fl '- a /(_l)2£/ 1 + lf/(_l) a Tl+«) > 


(4.33) 


where, in the second step, a' K l - a'j was used. After replacing (2a' K ^ + 1) by (2 a'j + 1) and properly adjusting the ordering of 
multi-products of V 2 a + 1 , we have 

2 ( -i)«/-i+f/+«A+A ( _i)«;-«; ( _i) 2 £ii j, x(ji,jjy-x(jj,j K y 
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1 ^7-1 Jj J 1 ^ w 7 

A~-l n-1 

[~[ V 2 a, + 1 P| V 2 a' s + 1 P| V 2 a t + 1 V27 + 1 


/-I 


1=2 


s=I 


k=K 


lJ lK 


. (4.34) 


The above result can be also directly written down from Eq. (4.24). In (4.24), the multi-products rim=y+i anc * summations 
2m=y +1 do not ex ' st for K = J + \. Hence the terms involve these multi-products and summations can be omitted, and this yields 
the result (4.34). 

The above discussions in case by case indicate that the general form of the operation can be written as 

~<JK, ,J> . " = V ( _ ] r ^ ]l+ a K ,j K{ _ , r -a' l{ _ , y, ( _ ,)-1£ +1 u X{ j^ j^Xijj, j K )k 

y * /A \ v / m\—mj--mj--m]c-m ll Z_J 

( a 'r a K- 1 ) 

(a/r jtc a K - i 
/ 7/J \ f Jk 


x -^(2a'j + l)( 2 a/ + 1 ) y j(2a' J + l)(2aj + 1 ) j ^ 1 J y 

y-i _ /. 

x [“[ y2a'+ l)(2a /+ 1 X-l) 0 !-' +a '- , + l <j J ' 

/=/+l ^ 

a:-i 

x [“[ V( 2<4 + l)( 2 a m + 1 )(-!)< 


+ 1 J Jl a i- J 

fl/ a/_l 


m=7+1 


+a,„_i + l J7m a m-l fl m 
1 a m a i n i 
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where the tuple a in the new intertwiner 


(tf J ) 
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i ■mi-mj- 
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M 







is given 

by 



r -mj-m K 

■••m„ 





• ,ai-u 

a'j, • • 


ax, ■ • ■ . 




(4.35) 


(4.36) 


For special cases of/ -I <7+1 and 77 - 1 < / + 1, the final results can be obtained from (4.35) by omitting the corresponding 
multi-products fl/j/'+i an d ]lm=/+i and summations Y/i=i+\ and 2 m=)+i- 

The second term in parenthesis of Eq. (4.2) can be calculated similarly to the former one. Here we omit the intermediate 
steps and directly write down the result as (a complete calculation is also shown in Appendix C) 


~<IJ;JK> 

‘Jijk 




M 


m\-"my-mj—mK"-m n 


Yj (- l) a, - 1+ ^ ,+aK+ j K (-l) ai ~ a, (— l)2b/u 7(_I) - jj)^X(jj, jx)^ 
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1 Q m flm-l 


j a 'j jj 

fly-i] 

ffly-i 

Jj 

fly] 

a a y-i 

jj ) 

l 1 

a 'j 

j J j 


X 




M 


(4.37) 


mymi—m 


Again, the final results for the special cases of / - 1 <7+1 and 77-1 < / + 1 can be obtained from (4.37) by omitting the 
corresponding multi-products n/="/+] and and summations Yjj=i+\ and Zm=y+i- 

Combining the results (4.35) with (4.37), for the case of 7 > 2 and 77 < n, the action of c/uk in Eq. (4.2) on the intertwiner 
can be explicitly written down as 

qiJK-ii? 3 ) M =V )«/-‘>'/C_i)I&l 1 A(_i)-Z^ + i7- XO - )1 
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/ ^ M 


jj a J-i\ \ a J-\ jj a J 


1 a', 


y-i 


Jj 


a j JJ 
(4.38) 


' mi —mi"-mj"-niK •••/«„ 


Again, we expect that the result (4.38) is general and also suitable for the remained special cases of 0 < 7 < 2 and K = n 
when «/ | and a K do not exist. While «/_i and a K do not exist in the intertwiner, we can ‘create’ them via (see Appendix B.3 for 
proof) 



(l) 2 ' 1 V2«o + 1 V2ai + 1 



Oi = j 1 


(4.39) 
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where dg = 0, d\ = 7l , and 7 is relabeled by a„. Then the intertwiner in (3.14) can be extended as 


(i v Jj ) M = (l) 2jl f] V 2 cn + 1 


/=0 


a 0 = 0 a 1= j i ' 


m n 

i* . (4.40) 

-+- 1 -»-M 

a n-l a n = J 


We are immediately awake to the fact that (4.38) is suitable for these special cases just by taking ao — 0, d\ - aj = jj, and 
= 7. 

Note that the operator quK at v only changes the intermediate angular momenta «/,••• , a E between I and K of the inter¬ 
twiner i v , but leaves the magnetic quantum numbers m i, • • • ,m n ,M invariant. The matrix elements of q IJK with respect to two 
given normalized spin network states T no ™(A) := X\ eeE(y \ -J2j e + 1 T :■ (A) read 
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where y is any graph bigger than y and y', T’fy) denotes the number of the edges in y, and d /Jn(g) is the Haar measure on S U(2). 
In the second step of Eq. (4.41), we notice that qu E at v only change the intertwiner i v , and if y differs from y', the integration 
with respect to the Haar measure gives a zero result. In the third step, the integration for the holonomies along the same edges 
but with different spins yields delta functions and the contractions of the intertwiner with its conjugate. In the forth and fifth 
steps, we have used the definition (A.33) of the inner product of the intertwiner space 7Yjj . associated to v. Thus the general 
expression (4.38) allows us to uniformly write the matrix elements of qjj E in the gauge-variant and gauge-invariant intertwiner, 
corresponding to resulting angular momentum J + 0 and 7 - 0 respectively, as 
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(4.42) 


where, in the third step, we have used the result (A.76). Note that the multi-products n/jz+i and !!«=/+1 anc l summations Y/i=i +1 
and 2^I ) +1 should be omitted for the cases of / < 7 + 2 and K < J + 2, and we need to set oq = 0, a 1 = a\ = j\ , and a„ = J when 
0 < 7 < 2 and K = n, respectively. Under exchange a <-» a', the expression in the square bracket of (4.42) is antisymmetric, 
while the other terms leave invariant because of the symmetric properties of 67 -symbol, symmetry of the delta function and the 
fact that (-1 ) a, ~ a ' = (-1 ) a i~ ai . Hence the matrix elements of quK are antisymmetric, i.e.. 


(a'\quK\d) - ~{a\qu K \a'). 


(4.43) 


The matrix element formula (4.42) derived in graphical method is the same as the formula obtained from algebraic manipulation 
for the case of 7 > 1 and / > 7 + 1 in [13], although different ways were adopted to deal with the recoupling problem. Moreover, 
as shown in above discussions, the formula (4.42) is also valid for other cases and hence can be regarded as a general expression. 
Finally, we consider some special cases which usually appear. With the following values of 67 -symbols [27], 
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where s = a + b + c, the general matrix element formula (4.42) can be simplified in the following special cases. 

(1)7= 1,7 = 2,77 = 3 

In this case, the general matrix element formula (4.42) reduces to 
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(4.47) 
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where we have used a {) = 0, a t = a\ = j\ for I = 1 in the first step, and (4.44) and (4.45) in the second step. Moreover, we can 
further simplify the result (4.47), since the triangular conditions on the 6 /-symbols will constrain the values of a' in (4.47) as 


fli = 


a 2 - 1 

a 2 

a 2 + 1 


(4.48) 


Denoting \a 2 ) = \a 2 ,a 2 , ■ ■ ■) and \a 2 - 1) = \a 2 - 1, a 2 , ■ ■ •), we get 
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where we have used the fact that (-1 ) 2jl + 2 i 2+2a2 (—] ) 2 / 3 + 2 a 3 + 2 a 2 _ ^ due to the triangle condition for (J \, j 2 , a 2 ) and ( 73, a 2 ,a 2 ). 

(II) / = 1,7 = 2, A - = 4 

In this case, the general matrix element formula (4.42) reduces to 
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(4.50) 


where we have used uq — 0, a\ = a\ = j\ for I — 1 in the first step, (4.44) and (4.45) in the second step, and the fact that 
a 2 + a ' 2 € N => (-l) 2 ( a 2 +fl 2 +1 ) = 1 in the last step. 

(Ill) 7 = 1,7 = 3 , A = 4 

In this case, the general matrix element formula (4.42) reduces to 
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(4.51) 


where ao — 0 and a\ — a\ = j\ have been used in the first step, and the 67 -symbol in Eq. (4.44) has been used in the second step. 
(IV) 1-2, J-3,K-4 

In this case, we have 0 /_] = ci\ — j\. Then the general matrix element formula (4.42) reduces to 
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(4.52) 

(4.53) 


5 Summary and discussion 

In the previous sections, the graphical method developed by Yutsis and Brink and their extensions, which suit the requirement 
of representing the holonomies and the intertwiners, are introduced to LQG. We firstly represent the algebraic formula by its 
corresponding graphical formula in an unique and unambiguous way. Then the matrix elements of the operator quK, which 
is the basic building block of the volume operator, are calculated via the simple rules of transforming graphs. Note that the 
calculations that we did by the graphical method can also be performed by conventional algebraic techniques. Also, to every 
graphical reduction, there is a corresponding algebraic reduction because of the correspondence between the graphical and 
algebraic formulae. However, it is obvious that our graphical method is more concise, intuitive and visual. 

Note that in our graphical representation, a gauge-invariant intertwiner associated to a vertex v of a standard graph at which 
n edges with spin 71 , • • • , j n incident is represented by 


n— 1 m, n 

[~[ \I-Qi + I [_-| 

/=2 

Taking account of Eqs. (A.43) and (A.50), the formula (5.1) is equal to (-l) 2 - 7 " times of the formula 


n—2 m ‘ m 2 m,,-! m„ 

n V 2 «; + 1 j \ j -\ 4 ■ 

1=2 ~ _ 



(5.1) 


(5.2) 


Since the only difference between (5.2) and (5.1) is a factor (-l) 2 *, the formula (5.2) is also used to represent the gauge-invariant 
intertwiner in literatures. 

Let us compare our calculation with those existing in literature. The operator qux can be represented by the following three 
forms [11] 


qim := -4^444 = 4 [4,4.44] - [(4 + 4/>(4 + A) 2 ] • 


(5.3) 


The first and the third forms (equalities) of the expression (5.3) were adopted as the starting points respectively in [12] and in 
[11, 13], and their matrix elements are calculated by graphical and algebraic methods respectively. In this paper, we considered 
the second expression (equality) of c/uk and derived its matrix elements by the graphical method introduced in section 3. In [12], 
to compute the closed formula, Pietri and Rovelli adopted the graphical Penrose binor calulus and tangle-theoretic recoupling 
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theory to deal with recoupling problems in the gauge-invariant context. However the rules of transforming graphs were not shown 
in [12]. So it is not obvious whether the rules uniquely correspond to the algebraic manipulation of formula. Note that the idea in 
[12] to employ the first equality of (5.3) to calculate the volume operator can also be carried out by our unique and unambiguous 
rule of graphical calculation. From (3.5), we have 


quK ■ 


U' d ) M = -\ie l]k J l ei J{/ eK • 

\ / m\—mi■••mj-m.fr-mn 1 K \ / m\—mi—mj—mK—m n 


- 


M 


mi •••m I "-mj---m K ---m n 


mi —mj—mj—m K —m n 


4 etjdnj.irdr m , [^(t,-)]"' mj [; n jK (r k )] m ' m , 
'(-40W Xj.W' m ,[^(Tv)]'^ m r^(T P )]< 


M 


m K ’ 


(5.4) 


where, in the last step, we have used the following identity (see Appendix B.4 for proof) 


e ijk [n h {Td] m ' mi [*j K (T k )] mK mK = -iWM^]"' 


m K ’ 


(5.5) 


with 


e pyp 



1 1 
y p, 


(5.6) 


Note that here one has e-i o+i = 1 (see also Appendix B.4 for proof). Notice that both [7^ ; ■(^ /J )]"' , and e /Jvp (given by a special 3 j- 
symbol) in Eq. (5.4) have corresponding graphical representations. The action of quK, corresponding to (5.4), on an intertwiner 
is given by 


m,i m.i +1 m.;_i m_, wij+i m,K -i m K 


n —1 

quK • J - [ yj2ai + 1 V2T + 1 

i=2 


— CLi — 0/+1 ' ’ ’ Oj 

V6 n ~ ] 

jj, j k) 2 ]~[ + 1 V2 J + 1 j, 

-2 — aj 1 

nr 

— Oj 

n/+i 

ji+i 

a J+ i • • • arc 

ij~i m 

jj- 1 

2 — Orr- 

j 

J 

nj+i m 

jj+i 

K — 1 T 

JK-1 

r 

lK 

JK 

+ 

’JK 

i= 2 

— flr+l • • • Oj -2 — 0.7-1 — 

oij — aj+i ■ ■ ■ o.k -2 — o/r-il 

1 J 



l l 


(5.7) 


where X(j,, jj , j K ) = 2j,(2j, + 1 )(2j, + 2)2jj(2jj + 1)(2 jj + 2)2j K (2j K + 1)(2 j K + 2). The derivation of the action of q IJK on the 
intertwiner in the graphical method is to remove the three curves with spin 1 in (5.7) by using the previous rules of transforming 
graphs. The identities in Eqs. (4.10), (4.11), (4.13) and (4.14) enable us to reduce the graphical formula (5.7) as 


y F(a,a! ,a" ,a r ") 


mi 

3 1 



(5.8) 


where the factor F(a, a', a", a"') involves the intermediate momenta a, a', a" and a'" in the intertwiner. By the following graphi¬ 
cal identity 



we can remove the three curves with spin 1 and obtain the final result, which coincides with (4.38). The closed formula of the 
volume operator was also derived by Brunnemann and Thiemann in [11, 13] using the algebraic techniques. The derivation 
process in [11, 13] is rigorous but rather abstract and awkward. Our graphical method is convenient and visual, and our result 
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(4.42) coincide with the formula derived by the algebraic calculation for the case of I > 1 and 7 > I + 1 in [13]. Moreover, our 
analysis shows that the formula (4.42) is also valid for other cases and hence can be regarded as a general expression. 

In principle, the matrix elements of qjjK in Eq. (4.42) enable us finally write down the action of the volume operator on the 
spin network states. We denote 


it 6 /I 3 

7v = ?(e/ ’ ej ’ eK) q,]K ' 

I<J<K, ejC\ejC\eK=v 

With the matrix elements of quK, we can get the eigenvalues and corresponding eigenstates of q v as 


(5.10) 


9vH*> ~ ^-q v \^-q v ') ■ 


(5.11) 


Then we can write down the action of V v on the intertwiner \i v ) associated to v as 


V v |i v > = VSj l*v» = 2 Vi§vi U«v)(^vl J 'v> = ^ <w 


I4>- 


(5.12) 


However, when the dimension of the intertwiner space associated to v is bigger than nine, we can not diagonalize q v analytically. 
This prevents us from explicitly writing down the whole formula for the action of V v . 


Acknowledgments 

The authors would like to thank Antonia Zipfel for helpful discussions. J. Y. would also like to thank Chopin Soo and Hoi-Lai Yu 
for useful discussions. J. Y. is supported in part by NSFC No. 11347006, by the Institute of Physics, Academia Sinica, Taiwan, 
and by the Natural Science Foundation of Guizhou University (No. 47 in 2013). Y. M. is supported in part by the NSFC (Grant 
Nos. 11235003 and 11475023) and the Research Fund for the Doctoral Program of Higher Education of China. 


Appendix A Elements of graphical representation and calculation 


A.l Representations of SU(2), Clebsch-Gordan decomposition, and the intertwiner 

To every non-negative integer or half-integer j (i.e., for j = 0, 1, |, • • •), there exists an irreducible representation nj of S U( 2), 

specified by j, on a Hilbert space 77, with dimension 2 j + 1. The orthonormal basis of 77, may be denoted by [e„ }, or {\jm)\ in 
Dirac’s notation, where m = - j , -j + 1, • • • , j. Given two irreducible representations ttj l and nj 2 of S U( 2) on 77, and 77 ; ,, the 
tensor product representation ttj [ ® n p on 77,, 877,, is (2 j\ + 1)(2 ji + 1 )-dimensional reducible representation of S U( 2). Clebsch- 
Gordan theorem tells us that the representation n J{ ® nj 2 can be decomposed into a direct sum of irreducible representations 7 xj 
on 77y, where J e {\j\ - j 2 1, • • • , j\ + ji }, formally. 


71+72 


1 J\J2 


{n h {g)®n h {g)){I hh r l = 0 tt j(g), VgeSU( 2 ), 


(A.l) 


J=\jl~j2\ 


where Ijj 2 is called the intertwining operator in representation theory of groups [33]. Given bases e\l\ ] ® e„ 2 6 77 ;i ® 77 j 2 (or 
\j]m ] j 2 m 2 ) = \jimi) ® \j 2 m 2 ) e 77,, ® 77y 2 ) and e { j^ e 77/ (or \JM) e 77y), the components of /„ j 2 , as matrix elements, are given 
by 




JM 


_ / (J) gO'i), 

' M ' I’h 


Si 2 ) 


\ M t < 

= e (jA e ‘ 


O i) o, £ Ui) 


\ / \ JM 

)• or Vhp) mm2 =(JM\j imi j 2 m 2 ), (A.2) 


where e^ J) is the dual basis or the basis of 77*. are called the complex conjugates of the Clebsch-Gordan coefficients 

(CGCs) or the intertwiners. In matrix form, the row and column indices of /, I /2 are denoted by the latter (upper) and former 

/ v. JM 

(down) indices JM and m\m 2 respectively. Our convention enable us to regard (//,/ 2 ) as components of the intertwiner 

tensor whose indices can be lowered and raised by a “metric” which will be introduced by (A.l). 

The representation n s of S U( 2) on 77, induces a conjugate representation n* of S U(2) on 77* via 




(A.3) 
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where [e'" ) is the orthogonal basis of 77*. Furthermore, the irreducible and unitary properties of n, are preserved to its conjugate 

\J) J J 

representation n*. If the representation ii j is unitary, there exists an unitary operator C (j] : 77, —> 77*, such that [28] 

C u >nj(g) = 7tj(gTC u > o C°^(g)C w_1 =7tj(gy. (A.4) 

The operator in fact defines an isomorphism between 77, and 77* by 

e m=C^ m ,e^, or \j m ) = C < ±,\jm')*, ( A.5) 

e ?i) = C U)' e m” or \jrn)* =C™’\jm'), (A. 6 ) 

where C™ = (C ur ) mn . The operator C ij) and its inverse 1 play an important role also in quantum field theories, whose 
components in the bases of 77, and 77* are given by [28] 


satisfying 


Obviously, c! and C™ ! satisfy 


cl := (-i y- n s n ,. m = (-1 ) j+m 6 m ,- n , 

C mn s (C <j)-' r (_1 y-’n 6m _ n = (-1)^"^ 

fiU) f^nm' _ (-im'nfiU) _ £in' 

/~i(j)/-inm' _ s-im'nsi(j) _ 

— ^(j) — v J-/ • 

c^, = (-i) 2 d, c™' = (-d 2 ^;-, 


(A.7) 
(A. 8 ) 

(A.9) 
(A. 10) 


(A.ll) 


which implies that C'l and C™" are symmetric for integer 7 , and anti-symmetric for half-odd integer j. The operator C <J) and its 

inverse C iJ> 1 can be used to lower and raise indices of the tensors on 77,. Hence Cl behaves like a metric tensor. Eq.(A.4) can 
be written in the form of its components as 


Cl, [Xj(g)] m ' n ,C% = [nj(gn 

The fact that the representation 7 r, is unitary implies 

xjig)' = *j(gr l = xjig- 1 )- 


(A. 12) 


(A.13) 


Hence 

ix,(g-')r m = i*j(gr l r m = [njign m = Wsyy\ = wi- n = ci m , , (a. m) 

where the overline denotes complex conjugation, and we have used (A. 3) and (A. 12) in the last two steps. By the map C (;) we 
can define a natural inner product on 77* as 

0 C U) f , C^g) Wj := (g, fy Hj , V/, g e 77,. (A. 15) 

Then the base transformation I J in 77,. ® 77 „ induces the base transformation /* . in 77* ® 77* by 

71J2 ' /1 J2 J 1J2 yi 72 J 

\ m i m 2 _ / \ JM _ m 2 m f 2 /j 

1 JuhJjM ~ V jljl )mum - 2-1 L 0'l) L 02> 

m'j 

For given 71,72 and 7, we denote (/ A j = which projects the bases ® of 77,, ® 77,-, onto of 

77y c 77,, ® 77 , 2 . The corresponding matrix elements of representations 717 , ® 7r, 2 and nj in the two bases, respectively, are related 
to each other by the so-called Clebsch-Gordan series 

Z W,AL« M[ v.(s)r„[v,(s)r-„((/j; A r , )„™ = (a.i?) 

mi,m2,ni,«2 


c u) 


M’M ’ 


(A.16) 
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Now let us consider the decomposition of the tensor product 7tj l ® • • • ® jtj n of irreducible representations of S U( 2) for n > 2. 
The composition involves n — 1 decompositions of the tensor products of two representations as (A.l) and the choice of the 
decomposition schemes. Denote rz, (i = 2, • • • , n) the irreducible representations that appeared in the i-th decomposition for a 
given scheme. In what follows, we consider the standard scheme where we firstly decompose jtj, <S> n j2 into (0) n tll , and then 
decompose n a , ® jtj 2 into Q) n a2 , and so on. We also denote a = (ci 2 , ■ ■ ■ For given j\ ■ • • , j n , allowable J, and compatible 
vector a = ( 02 , • ■ • the corresponding Clebsch-Gordan series read 




m !,••• ,m n ,n ],• 


(A.18) 


where 



m.\—m n 
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(/? . ) 

V Jv-JnJn 


JM 
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/r \ ^ . 

/ \ fli+1^1+1 

It \ 

JM 

\ •^ 2 /m 1 m 2 



r . n -\m n 


are the general CGCs, the intertwiners, and often rewritten in quantum mechanics in the form 

(JM;d\jimij 2 m 2 ■ ■ ■ j„m „> = ^ {a 2 k 2 \j\m { j 1 m 2 ){a 2 k- i \a 2 k 2 j 2 m ? ,) ■ ■ ■ {.JM\a„-ik n -ij n m „>. 

k 2 ,-,k n . 1 


(A.19) 


(A.20) 


It is easy to generalize the above results to the decomposition of the tensor product nj 1 ® • • • ®7tj k (g) ® jtj k+1 (g) _I ® • • • ® 7H„(g) _1 of 
k repsentations and /; - k inverse of representations into a direct sum of irreducible representations ttj on 'Hj. The corresponding 
Clebsch-Gordan series read 


2 ( 7 ,W X 








ll«A:+l 


wAIO'), 


m 1,— ,m„,» 1,— ,n n 


[X /(*)]%. 


(A.21) 


where 


> Jy"Jn/m\—mk 


tjJ-,3 \ M ^-.nt+imt+1 . . . Qn„m„ 

V ji"'jn/mi—mkni [ +i—n n (jk+i) (jn) 


(A.22) 


Eq. (A.21) can be written as 


E w'r'wT, 


mk+\--m n M 


[x h (g)r ni • • • = 2 bcte)]" 


n k +\—n n N 


(A.23) 


which, in the case of / = 0, reduces to 

E [”».>WT- 1 


>nk+r-m„0 


i=Ky,.' a - 24 > 


Hence the tensor 
introduce 


(/? : f ) 

V Jl-JnJn 


m£ + i"-m„0 


is also called the invariant tensor. For the convenience of graphical representation, we 


M := ct\hnujimz ■ ■ • ./>»„> • (A.25) 

\ 71 Jn/m\—m„ \ 71 Jn/m\—m„ 


Notice that the factor (-I ) 71 7i 7 in Eq. (A.25) involves only the spins j\, - ■ ■ , j„ and J, not the intermediate momenta 
«!,••• ,a„-\. From now on, the inertwiners refer in particular to z.^E . The Clebsch-Gordan series (A.18) can be written in terms 

of i J ’ “. as 

h-Jn 




mi,— ,m„ ,«!,••• ,n n 


and its inverse reads 


- E (0'%*”. 


J.M,N 


(A.26) 


(A.27) 
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In the special case of n = 2, the Clebsch-Gordan series read 


[njM)r ni [n j2 (g)r-„ 2 = 2 


(A.28) 


J,M,N 


The fact that the operator i d . is unitary and its matrix elements take real numbers results in 


- (or - msjl 


(A.29) 


Similarly, the relation (A. 16) can be generalized as 


K'/I Iw - ('/ /,“ Zj c 0i) ”' C Un) \ l h-in) m \-m' C M'M - \ l jv’ (A.30) 


which shows that the GCG is real (see Eq. (A.75) for proof with graphical method). Moreover, the Clebsch-Gordan series (A.21) 
for the general case can also be written in terms of i a ’. as 


y (i j;a ) 

2_J V JV-JnJmy 


m k +\—m n M 


mi,— ,m n ,ni,— ,n„ 


1^(8)]"" „, • • • [Xj k m m \ lk [n h MT l T M mt+1 • • • [*/„(*) Ov,....,,/' 1 

(A.31) 


where 


(i**. ) 

V JvJnJn 


lowing inner product 


rrik+\—m n M / 

i J ’ a ) 


\-m k V 

ivin) 

m\-m k n k+ v 

, the intertwiner 

space 

Wju-J, 

((i J '’ a ') M ' . 

(i J ’ a 

X M' 

y\ J\"']n)m\•■■m n 

' yji-j. 
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U (A + l) ' ' ' ^0„) 


(A.32) 

fif’".) with the fol- 

\ Jv'Jntm\-m n 


M 


Ji’-’jn mi,—,m n 


= y M '(i j ’ a ) 1 

Zj 1 ./ /• m Vti-tJ, 


mi,— ,m n 

= 5jj’5M,M’5a,a' 


(A.33) 


( 7 . ji \ M 

ij’ j j is unitary. The 

second one will be shown in Eq. (A.76). 


A.2 The basic components of the graphical representation and simple rules of transforming graphs 

In this subsection, we introduce the basic components of the graphical representation and simple rules of transforming graphs 
[20], A graphical representation for the matrix elements of irreducible representations of S U( 2) is also proposed. A graphical 
representation is a correspondence between graphical and algebraic formulae. Each term in an algebraic formula is represented 
by a component of an appropriate graph in a unique and unambiguous way. 

The Wigner 3 /-symbol is associated with the coupling of three angular momenta to give zero resultant. The 3 /-symbol has 
simple symmetric properties and hence is easier to be handled than the CGC. The 3/-symbol is defined in terms of the CGC by 
[34, 28, 27] 


m, m 2 m] ~ (-iy W 2 +; 3 (-l ) h+m K2h + l ) 4 {jim.hmilh - m 3 >. (A.34) 

The 3/-symbol takes non-vanishing value when the parameters of the upper row ( j\,j 2 , 73 ) satisfy the triangular condition (i.e., 
|/i — / 2 I < / 3 < / 1 + 72 ) and when the sum of the parameters of the lower row (ni \, m 2 , m 2 ) is zero. The parameters /; and m, are 
simultaneously integers or half-integers, such that each of the following numbers 

ji + mt, ji ~ mi, j\ + /2 + 73 , -j 1 + 72 + . 73 , j\ - h + h, j\ + h - h (A.35) 

takes some integer. 
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The graphical representation of the 3 /-symbol was collected by Yutsis in [29] and slightly modified by Brink in [20] for 
convenience. The 3 /-symbol is represented by an oriented node with three lines, which stand for three coupling angular momenta 
71 , 72,73 incident at the node [20]. The orientation of the node is meant the cyclic order of the lines. A clockwise orientation is 
denoted by a ” sign and an anti-clockwise orientation by a “+” sign. Rotation of the diagram does not change the cyclic order 
of lines, and the angles between two lines as well as their lengths at a node have no significance. The 3 /-symbol can be written 
in the graphical form 


/ 7i 72 73 

l mi m 2 m2 




(A.36) 


The 3y'-symbol has the following properties. An even permutaion of the columns leaves the numerical value unchanged, while 
an odd permutation is equivalent to a multiplication by (-iy i+ ^ 2+j3 , i.e., 


p' 

72 

73 \ 

— r_]_y’i + 72+73 1 
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\m 1 

m 2 
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m 3 

m 2 / 



(A.37) 


Moreover, the 3y'-symbol has the symmetric property 
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_ (_] y/i+ 72+73 I 7l 
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\mi 
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-m 3 ) 


The orthogonality relation for 3y'-symbol is expressed as 
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\m 1 
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m 3 ) 




Furthermore, the 3y'-symbol is normalized as 
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(71 

72 

73 \ 

f7' 

72 

73 \ 

\W| 

m 2 

m 3 / 

\W| 

m 2 

m) 


= 1. 


(A.38) 


(A.39) 


(A.40) 


The “metric” tensor in Eq. (A. 7), which occurs in the contraction of two 3y'-symbols with the same y values, is denoted by 
a line with an arrow on it as 


C (j) =(-l V~ m S >=(-lV +m 'S > = 

m'm V L ' u m,-m V 1 ) u m,-m 


(A.41) 


and its inverse in Eq. (A. 8) can be expressed as 


C n,n' = (_1)7-^, = { -iy+ m ' 6i 


m-m' ~ m - 


The special 3y-symbol with one zero-valued angular momentum is related to the “metric” tensor by 


7 / O' 
m m! 0 
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A line with no arrow represents the expression 


V^./ +1 


°m,m' ~ m - 


The property in Eq. (A. 38) means 
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(A.42) 


(A.43) 


(A.44) 


j I r 1 2 z 3 - 

' I C (/l) L (72) C 0'3) 
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72 73 \ 

m 2 W 3/ 



/ii' ,m' ,mt 


_m 2 m; _m 3 m' 
C 0'2) C 0'3) 



72 73 \ 

m 2 m 3/ ' 

(A.45) 


In graphical representation, two lines representing the same angular momentum can be joined. Summation over a magnetic 
quantum number m is graphically represented by joining the free ends of the corresponding lines. Eq. (A.34) implies that the 
CGC can be represented graphically by 


m, m 2 


<7i«M 72'«2[73'»3> = (-l) jl 12 J3 V 2 73 + 1 


= ( _ 1 yi-72-73(2 7 - 3 + 1 ) i/ 2 ^ 


(A.46) 


Therefore, the graphical representation of the intertwiner defined in Eq. (A.25) is 

n -1 


m, m 2 


teii, : = n _ t 


i =2 

and the generalized intertwiner in Eq. (A.32) can be repressed by 

n— 1 


m 2 m k m k _ 


V Jy-JnJn 


m k+ \-m n M 


:= ]~[ V2fl/ + 1 V27 + 1 


3 2 

jk 

Jk+l 

' jn 

v 





■ Af 


(A.47) 


(A.48) 


As an ingredient of the flux operator, the spherical tensor [nj(T M )] m ' m in Eq. (3.7) can be represented graphically by 




V 2 ;( 2 ;+l)( 2 ; + 2 ) 

A m!- 
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7 + J 
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l -^2j{2j+l)(2j + 2) 



(A.49) 


Now we outline some useful rules of transforming graphs, which can simplify our calculation of the action of operators on 
the quantum states in LQG. The frequently used rules for adding or removing arrows in a graph read 


TO- 

TO- 

TO- 




(A.50) 
(A.51) 
(A.52) 

(A.53) 


which correspond to the algebraic formulae in Eqs. (A.9), (A. 10), (A.l 1) and (A.45). The rule to remove a closed loop in a graph 
reads 


n 

m 3 —— h 


7q / 3 n~ 

— = - m 3 -— 


273 + 1 


(A.54) 


which is related to the orthogonality relation for 3 /-symbol in Eq. (A.39). 

Coupling four angular momenta to a zero resultant will involve the /m-coefhcients. The /'m-coefficients corresponding to 
different coupling schemes are related by the 6 /-symbol. The 67 -symbol is defined by ([27] p. 94) 


7i 72 73 

74 75 76 


Z 

all m,m‘ 


7l 72 73 \ / 7l 75 76 1 / 74 72 76 W 74 75 73 

\m\ m 2 tn .5 m' 6 )\in' 4 m' 2 m' 5 m' 3 J O'O O 2 ) (h) (M) 0's) O'e) 


(A.55) 

The four triangular conditions satisfied by the six angular momenta in the 67 -symbol may be illustrated in the following way 

(A.56) 
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The 6/-symbol has the following symmetric properties. It is left invariant by any permutation of columns and also by an inter¬ 
change of the upper and lower arguments in each of any two columns, e.g.. 


7 i 7 2 7 3 
74 75 76 


Graphically, we can express the 6/-symbol in Eq. (A.55) as 
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(A.57) 


7i 72 73 

74 75 76 


76 72 74 

73 75 7i 




(A.58) 


where in the last step we have used Eq. (A.53) to remove three arrows. Taking account of the definition of 6 /-symbol in Eq. 
(A.55) and the fact that the 3y-symbol is normalized in Eq. (A.40), we can easily show the following identity (see [27] p.95) 
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Hence we have the following graphical identity 
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(A.59) 
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m 2 = ji 



(A.60) 


; schemes 
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72 

76)/ 

1/4 

73 

/5j\ 


76 \ c '«(K / 76 73 74 

F72f, / Os) m 3 m 4 


(Aril) 


corresponds to the following rule of transforming graphs 

mi 7714 

\Ji u/ „ 

= V (2 y 6 + 1)(- 

/J3 j2\ ; 

m 3 m 2 76 

Using Eq. (A.37) and the symmetric properties of 6y-symbol in Eq. (A.57), we can get that from Eq. (A.62) 


_l)7 2 +7 3+ 75+76 J 7l 72 76 

1/4 73 75 


mi 

Vi 74/ 

/O 2 73 s 

m 2 


48 £ / = 


m 3 


= V(2/6 + 1 )(-l J22+73-/5-71 jri 72 761 

V U4 73 75J 


m 2 76 


77t4 mi 

V 4 71/ 

. . . r 

74 73 75 J A 

m 3 m 2 


Similarly, using Eqs. (A.37) and (A.63), we have 


mi . 7714 

4>—= £(2 J6 + 1)(—1 )7'-7 2 +73+./ 4 |7i 74 76 


m 3 


m 2 76 


72 73 75 


m 3 . . mi 

\73 7i/ ‘ 


m 2 ‘ 


(A.62) 


(A.63) 


(A.64) 


7714 


Note that Brink’s graphical representation in [20] does not involve how to represent graphically the matrix element of the 
representation of S U{2). Here, we will extend the Brink’s representation and propose a graphical representation for the unitary 
irreducible representation nj of S 7/(2). The matrix element [Xj(g)] m n is denoted by a blue line with a hollow arrow (triangle) in 

Ms)r„=: m _ ^ (A.65) 


it as 
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The orientation of the arrow is from its row index m to its column index n. The matrix element [nj(g 1 )]" in Eq. (A. 14) can be 
presented by 

- n -- m = n—< - fA.66) 

Up to now, we have expressed the quantum states (the spin network states in (2.8)), and the two elementary operators (the 
holonomy and flux operators in (2.12) and (2.13) ) of LQG in the graphical form. The Clebsch-Gordan series in (A.28) can be 
represented by 


m i 



m 2 



m[ 

m' 2 


Yi 2 h + 1) mi 


- {> 

Pa 
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m' 2 


<=Y(2J3+V) mi 

h 


. . h 

^ - |> 


h-h 

32 
m' 2 


(A.67) 


where, in the second step, we used (A.52) to flip the orientations of two arrows. Using Eqs. (A.66) and (A.67), we have 
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= X/ 2 ' 3 +!) 

mi _ 

-fr- 

2 ' 3 -ii , 

- *“1- m l 


(AM) 
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m 2 
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A.3 The spin network states and their orthogonality in graphical method 

In this subsection, we will give a graphical representation for the spin network states and show that these states are orthonormal 
to each other by the graphical method developed in the previous two subsections. The spin network state consists of a standard 
graph y, the non-trival spins j associated to all edges of the graph, and intertwiners i to all vertices of the graph 4 . The normalized 
spin network state takes the following explicit formula 

rpA) := (g) i v • (g) 4 2 J e + 1 n h (h e (A)) ■ (g) , (A.69) 

veV(y) eeE(y) ve V(y) 


where ■ stands for contracting the upper (or former) indices of representation matrices nj e (h e (A)) with the indices of intertwines 
i v at true vertices v, the lower (or later) indices of nj e (h e (A)) with indices of the conjugate intertwiners i*- at pseudo vertices v. 
Comparing to the form (2.8), the state (A.69) is normalized. The scalar product of the spin network states is defined by the 
formula 

Lnorm^nomj ;= f FT du H (h e )T™?^T n0 ™AA), (A.70) 

V 7J.I /-Wkjn Js£/(2)W7)I^ £( 4) ^ YJ ' y 


where y is any graph bigger than y and y', |E(y)| denotes the number of the edges in y, and d j3n(g) is the Haar measure on S U(2). 
If y differs from y', e.g., the edge e' with spin j' e , belongs to y' but does not exist in y, then the orthogonal relation 


f d fi H (g) [T:j(g)rMr(8)f n ' = , 

JSU(2) Z J + 1 


tells us that the corresponding integration in (A.70) becomes 

f dfiH(he’) [ny,(he ')]'"', = 0 . 

JSU( 2) ' n 

Hence the nontrivial result corresponds to the case y = y'. Thus (A.70) reduces to 

/^norm ^norm \ = § f Tf A n H (h e ) 7™™(A) 7™™ (A) 

\ y,j,i /'Han JSU( 2)W»l e ^) ^ ^ 


(A.71) 


(A.72) 


4 The relation between the spin network state on the standard graph and the corresponding one on the original graph is discussed in section 3.2 
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(A.73) 


= Sy,y> f FI d/j H (h e )T™™(A)T no ™(A), 

Jsua^w 


Jsuvrw e t E i r) 

where, in the second step, we have used the fact that the Haar measure is normalized. By Integrating over all representation 
functions on the edges, one can obtain the contract of the complex conjugation of intertwiners with the corresponding intertwiners 
at vertices. Thus we have 


/ r norm r norm \ = 6 y y FI 6j f FI V ) "(f 7 ""') ", 

\ 7 ,jj Y J' J' 1 I 1 1 / j V JV"Jn)m v -m n V Ji"'Jn)mym n 


(A.74) 


eeE(y) 


veV{y) mi,-- ,m n 


where n is the number of the edges incident at v, and the complex conjugation of intertwiner related to the corresponding 
intertwiner (A.30) satisfies 


V Jl-Jn/n 


M 


Z r mim[ r m n m' n / . J- 3 \ M r (J) 

L Ul) ' ' ' L 0„) \ l h-in)m\ m' n M’A 

fflp- 

n — 1 7* 

| ] V 2 ^/ + 1 V2/ +14 

i =2 
n— 1 

= Y\ V 2a / + 1 V2/ + 1 

1=2 

= 0 **.) M , 


(A.75) 


where, in the third step, we have used Eq. (A.53) to remove arrows. Equation (A.75) reflects the fact that the CGs (and thus 
GCGs) are real. The sum in Eq. (A.74) can be calculated by graphical method as 


y (;**.) v-) _ y y-) y 

V \ Jyin)my--m n V JvJn)m\--m n V JvJn/i 

m\,—,m n m\ 

«-l ( _ 

= n ^(2a, + l)(2a; + 1) V(27 + 1)(2/' + 1) 


M' 


<-l + J' 


1=2 


n-1 



- [y 5a i ,a' i 6j,r5 M J 

i =2 


(A.76) 


where, in the third step, we have used Eqs. (A.54) and (A.50). Hence, the inner product of the spin network functions reads 


/ \ " _1 
/ 'T’norm 'r-norm \ _ c 1 | c 1 | 1 | c c c 

r yj7 ’ T y ,]'? 1 1 I 1 1 1 6 a^Sjj'6 M , 


eaE(y) 

= 1 


veV(y) i=2 


(A.77) 


Note that the orthonormal result we obtained in Eq. (A. 77) depends on the premise that the intertwiners at the same vertex v 
involve the same coupling scheme. If different coupling schemes at the same vertex were chosen, certain additional multiplication 
of 6/-symbols would appear in the result. If the spin network functions are gauge invariant, corresponding to J - 0 and M - 0, 
the two Kronecker delta functions become 5jj' = 6m,w = 1, which will not appear in the expression (A. 77). 

Appendix B Proofs for some identities 

B.l Proofs of algebraic identities in Eqs. ( 3 . 7 ) and ( 4 . 5 ) 

By the definition of the matrix element 


Mr,)]" 


f=0 


[nj(e tT ')V 


(B.l) 


32 
































we get 


M T l)] m m = ~2 V)0' + 1) - m'(m' + 1) (5rn',m-l - 2 ViO' + 1) - '"'(ill' “ 1) <W,m+l , 
[7r/T 2 )]'"' m = ^ V)(; + 1) - in'(in' + 1) (W,_| - ^ ViO' + 1) - - 1) <W,m+l , 

[^■(T 3 )] m ' m = —i III 6 , n 'j n . 

Hence the matrix elements of r ;( (/./ = 0, +1) defined in Eq. (3.6) read 
[7r/T 0 )r m = In Xr 3 )] m m = 

• V2 


Wr * ,)r - = T (Wr ‘ )r " +iWr2)r '-> = +i ^ ' , ^ +1) “ m ' (m '- 1) 4 *-~'' 

= +d=f[3Tj(ri)] ffl m -i[^j( T 2)r ra ) = -i-J- Vj(J + 1) - m'(m' + 1) o, l; „ . 


Taking account of the specialized formulae for the CGCs, 


/ m l —:- 

(jml0\jm > = . 6,„', m - - V 2 ./ + 1 


j 1 j 

m 0 m' 


1) 

/■ in- /\ I JO + !) - “ !) x h - , ( 

\ 2 j(j +1) V 

z ■ , , x , I JO + 1) - m'(m' + 1) j l j 1 j 

0ml - 1 |jrn ) = +J - - — - = ~ 0 2 J + 1 I. 


j l j I /■■>>!' m' 

m 1 m") O') 


£<m m 

m — 1 m" I 0) 


we get 


M^)] m ' m = l ~ y/2j(2j+l)(2j + 2) 


j 1 j 

in jj in' 


C^”' = l ~ y/2j(2j + 1)(2 j + 2) 


1 j j 

H in" m 


C 
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0) 


By definition (3.6), we have 


Then 


1 i 

T 1 =-pCr+i-T-!), T 2 = —(T +1 -t-T-O, T 3 =T 0 . 

V2 V2 


TjTj = Tin + T 2 T 2 + T 3 T 3 = —T+lT-l - T_lT+l + To To . 


Thus 


(B.2) 

(B.3) 

(B.4) 

(B.5) 

(B.6) 

(B.7) 

(B.8) 

(B.9) 

(B.10) 

(B.l 1) 

(B.12) 

(B.13) 


[xjMdr JnjMdr mj = -bH,( t+ 1 )] m ' m ,[njj (t_ 1 )]'”•' my - [^(T-Or^l^i)]”^ + [^(to)]"^^^)]”^ 


B.2 Proofs of graphical identities in Eqs. ( 4 . 10 ), ( 4 . 1 1 ), ( 4 . 13 ) and ( 4 . 14 ) 

Eq. (4.10) can be proved by 


(B.14) 
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(B.15) 
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where Eq. (A.37) was used in the first and third steps, (A.62) was used in the second step, and we used the fact (—= 1 and 
the symmetric properties (A.57) of the 6y-symbol in the last step. 

Eq. (4.11) can be shown by 
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dK-1 

k/T-i 

1 

dK 


(B.16) 


where Eqs. (A.62) and (A.52) were used in the second and the third steps, and in the last step we used the fact (-1) = 1 and 

the symmetric properties (A.57) of the 6/-symbol. 

Eq. (4.13) can be proved by 
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1 O-l/ \f L l -1 a l 
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(B.17) 


where Eq. (A.37) was used in the second and fourth steps, (A.62) was used in the third step, and in the last step we used the 
symmetric properties (A.57) of the 6/-symbol and the exponents were simplified. 

Eq. (4.14) can be proved by 

V.-i 1 / ,, , \a„—i b' ri 

— = V-fc-4 = (-i) am+fc ” +1 

b "> Am 
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fr - * l 7m a m I 
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= y(2^_j + i)(-i) fc :- 1 + ^-i +1 (-i)^-^ 
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// . 
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d m 

*2/77-1 


^TO-1 — 
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(B.18) 


where Eq. (A.37) was used in the second and fourth steps, (A.62) was used in the third step, and in the last step we used the 
symmetric properties (A.57) of the 6/-symbol and the exponents were simplified. 
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B.3 Proofs of graphical identities in Eqs. ( 4 . 18 ) and ( 4 . 39 ) 

Eq. (4.18) can be proved by 


a!j ~ .“J ~ „ . , . , a'j + aj - b'j „ 

K J -^ ^ J—► .1 — ^+ l Kj -^ * J— *— h _ ^ Y y 


cij+a'j+1 k 'j- 




~l'j _ ^_j^ay+fly+l kj- 
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■ Ki » 


= (_l)'u +a / +1 (_l) 2a / k 'j —^“r— j —— 1 'j - (_i y'j- a J+i — a '' bj k , _ 

2 a'j + 1 J 

where identities in Eqs. (A. 37), (A.53), (A.50), (A.51) and (A.54) were used from the first to last steps. 
Eq. (4.39) can be proved by 


(B.19) 



y/2ji + 1 
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V^/i + 1 J j = y/2ji + 1 

o-» _!_ « . - 
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0 _ 3i 


= (l) 2jI V2;'i + 1 * = d) lh y/2a 0 + 1 V2«i + 1 A , (B.20) 

o ►— 1 —► o ►t-4— *r ~-— 

0 Ji a 0 = 0 a, = j! 

where the rules (A.50), (A.43), (A.53), (A.50) and (A.52) were used from the first to fifth steps, and in last step we denoted 
ciq = 0, a\ = j\. 


B.4 Proofs of algebraic identities in Eqs. ( 5 . 5 ) and ( 5 . 6 ) 

To prove Eq. (5.5), we denote [nj I (T i )] ni mi [7t jj (Tj)] Hj mj [7ij K (T k )] ni mK by [t.-T/T *] n , m , njmjriKmK . Then we have 

eijk[TiTjT k ] mminjmjn K mK ~ ^1231 2 "I '^2^3\n i minjmjn K m K + ^l'i2\j\^^l\niminjmjtiKmK ^213\T2^1^3\niminjmjn K mK + Ql2[ T 3 1 "l T 2]n;m/)ijmj«s:ms- 

"t" ^231 [t"2T3T] jn^m p ^321 [t*3t"2f I ]n/m/njmjn p tUF; 

— [r \T2 ~ 2" 2'^\\nimjnjmj\T'i\nxmK [Tsl/i/m/lTlTo — TjT | \n/nijn p m p ~ [T | T 3 T 2 T2T^T\\n I m I njmjn^mf' - IB.2 1 ) 

Taking account of (B.12), we have 

1 i 

t it 2 - r 2 ri =- — — [(r+i - t_i)(t +1 + t_i) - (t +1 + t_i)(t+i - t i)] = -/(t+iT-! - r_iT + i). (B.22) 

V2 V2 

Hence, we get 


e i j k [n h {TdT I m^ j AT j )r mj [n jK {T k )T K m K 

— [Tl^"2 T2'T\\niminjmj\T'i\nKmK [ 2~3 ] /?///; / [t"l To — T 2 T \ ^njmjn p m p ~ |r | T 3 T 2 — T 2 T 3 T [ ^nim.!njmjn p m p 

= - i[r+iT 1 - T-iT +1 ] niminjmj [Toin K m K ~ i[To\ ni m, [T+lT-i - T_iT +l] nj mjn K m K ~ (“OlT + 1T 0 T_1 T—iToT+i^nimjfijmjtiKmK 
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= _ t|A+l-10 T +l T -l T 0 + f-l + 10 T -lT+l T 0 + 6 o + l-lTo T +l' r -l + e0-l + l T 0 T -lf+l + f-l 0 + 1' r -l T ()T+l + e + io_iT + iToT_i]„, m;njmjns . mj . 


— pT y Tp\mminjmjnKinK 

where e pvp is the Levi-Civita symbol defined by e- 10+1 = 1. 

To show Eq. (5.6), notice that 

1 1 A J_ 

-i o ij- vr 

Then one has 

1 1 l\ 1 

-1 0 lj“V6 e_1 ° +1 


(B.23) 


(B.24) 


(B.25) 


Recalling the symmetric property of the 3 /-symbol in Appendix A.2 , an even permutation of the columns leaves the numerical 
value unchanged, while an odd permutation will lead to a factor (-1) 1+1+1 = -1 for the 3/-symbol in Eq. (A. 37). These 

symmetries of the 3/-symbol are the same as those of e pvp . Hence we have e jJvp = V6 


1 1 1 

P v P. 
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Appendix C Complete calculations for the action of the second term in the paren¬ 
thesis of Eq. (4.2) 

The action of the second term in parenthesis of Eq. (4.2) can be explicitly calculated as 

M 




n -1 

: X(ji, jj) 2 X(jj, j K ) 2 J - [ y/2a, + 1 V2/ + 1 


raj-i m j m,j+i m K _i mu 
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(C.l) 


ai — a I+ i ■ ■ ■ aj _ 2 — a.;_i| — I aj — aj+i ■ • ■ ok - 2 — clk- 


Similar to the calculations in section 4 for the first term in the parenthesis of (4.2), we perform the calculation of (C.l) in four 
steps. In particular, from the second step we first consider the case of J > I + 1 and K > J + 1. In the first and second steps, the 
calculation are completely similar to those for the first term. Thus we can write down the result directly as 
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(C.2) 


In the third step, after dragging two nodes attached to the jj line down to the horizontal lines by using firstly (4.10) and then 
(4.1 1), we have 
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Now we turn to the forth (last) step. By using identity (4.18) and summing over b' . and a',, we obtain 


Yj X + lx-iy 1 '- 1- ^ | ^ ^ a ! \ x (2/7^! + i)(_i)^-^-i + ^ +i J a * j K aK ~ x 
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The exponents in (C.5) can be simplified as 
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By (4.21) we obtain the result 
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In what follows we now consider the remaining cases. 

(I). The first case of / = I + 1 and K — J + 1 

The result can be directly obtained from that in the case of / > I + 1 and K > J + 1 in Eq. (C.5) by omitting those terms appeared 
in second step as 
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The exponents of (C.8) can be simplified as 
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(C.9) 


where the identities aij \ = a/ and a ^= a'j were used in the second step. After replacing (2a’ Kl + 1) by (2 a'j + 1) and properly 
adjusting the ordering of multi-products of V2a + 1, Eq. (C.8) reduces to 
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(C.10) 


In fact, (C.10) can also be regarded as a special case of Eq. (C.7). In (C.7), the multi-products n/=y+i anc * rim=y+i ar| d summa¬ 
tions Y/i=i+ i anc l 2m=y+i do not exist for the case of / = I + 1 and K - J + 1. Hence the terms involving these multi-products 
and summations can be omitted, and this yields the result (C.10). 

(II). The second case of / = I + 1 and K > J + 1 

The result can be directly obtained from that in the case of J > I + 1 and K > J + 1 by omitting certain terms appeared in the first 
line of the second step in Eq. (C.5) as 
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The exponents in (C.l 1) can be simplified as 
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(C.12) 


where in the second step we replaced aj-\ by a /. By using (4.21) and properly adjusting the ordering of multi-products of 
V2^TT , we have 
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(C.13) 


Again, (C.13) can also be directly written down from Eq. (C.7). In (C.l), the multi-products nf=/+i an d summations 2/4 /+1 do 
not exist for the case of / = I + 1. Hence the terms involving these multi-products and summations can be omitted, and this 
yields the result (C.13). 


(III). The third case of / > I + 1 and K = J + 1 

The result can be directly obtained from that in the case of / > I + 1 and K > J + 1 by omitting some terms appeared in the 
second line of the second step in Eq. (C.5) as 
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The exponents in (C.14) can be simplified as 
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(C.15) 


where in the second step we used a’ K -\ ~ a 'r By replacing (2 a' K _ x + 1) by (2a) + 1) and properly adjusting the ordering of 
multi-products of V2a + 1, we have 
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The above result can also be directly written down from Eq. (C.7). In (C.7), the multi-products n,„=j+i and summations YJm=j+ i 
do not exist for the case of K = J + 1. Hence the terms involving these multi-products and summations can be omitted, and this 
yields the result (C.16). 

The above discussions in case by case indicate that the general form of the action can be written as (C.7), i.e., 
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For special cases of J— 1 <7+1 and K— 1 <7+1, the results can be obtained from the form (C. 17) by omitting the corresponding 
multi-products fl/wAi and Wfn^j+\ and summations Yj=i+i and Em=)+ 1 - 
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